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ABSTRACT 


RUSAP is a FORTRAN IV computer program designed to determine the 
performance characteristics (power-to -weight ratio, blanket tension, struc- 
tural member section dimensions, and resonant frequencies) of large-area, 
roll-up solar arrays of the single-boom, tensioned-substrate design. The 
program includes the determination of the size and weight of the base struc- 
ture supporting the boom and blanket and the determination of the blanket 
tension and deployable boom stiffness needed to achieve the minimum-weight 
design for a specified frequency for the first mode of vibration. This report 
provides a complete listing of the program, a description of the theoretical 
background, and all information necessary to use the program. 
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I. INTRODUCTION 


In recent years, considerable emphasis has been placed on the devel- 
opment of large-area solar arrays with high power-to-weight ratios and 
small packaging volumes. One of the concepts currently being developed is 
the single-boom roll-up array shown in Fig. 1 and described in Ref. 1. 

The design consists of two flexible cell blankets tensioned between 
spacecraft-mounted storage drums and a leading-edge beam. The array is 
erected by a deployable boom that is connected between the leading -edge 
beam and its supporting structure on the spacecraft. 

In studying potential applications for this design, one must predict the 
performance characteristics of arrays of widely varying sizes and natural 
frequencies. Because of the complex nature of the relationship between the 
size of the array, its first-mode natural frequency, and its structural 
parameters, it is difficult to predict the performance characteristics of 
arrays significantly different from the engineering prototype. To help solve 
this problem, a computer routine RUSAP was programmed to calculate the 
optimum size for structural members and the optimum blanket tension for an 
array that is to have a given area and width and a specified lowest deployed 
natural frequency. The program is composed of two basic parts: (1) an 

analysis of the deployable boom and array blanket, and (2) an analysis and 
sizing of the base structure supporting the boom and array blankets. 

II. THEORETICAL BACKGROUND 

A. Optimum Boom and Blanket Tension Analysis 

If the power, area, and width of the array are specified, the length 
(and thus the size) of the array is fixed. Because the weight of the cell 
blanket and support structure is essentially fixed for an array of a given 
size, the relationship between total weight and deployed natural frequency is 
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almost entirely controlled by the weight of the deployable boom. This 
implies that, when the size of the array is specified, the optimum array for 
a specific lowest deployed natural frequency can be defined as that with the 
lightest boom. 

Because there are two possible fir st-vibration modes for a deployed 
array (Fig. 2), the cross-sectional size (and thus the weight) of the boom is 
determined by the following requirement: the boom must have sufficient 
stiffness to maintain the first-bending frequency equal to or greater than the 
required minimum frequency while loaded by sufficient blanket tension to 
maintain the first-torsion frequency at or above the minimum. Figure 3 
shows the variation in first-bending and first-torsion frequencies as a func- 
tion of the blanket tension for a single-boom array. As is shown in this fig- 
ure, the optimum tension for a particular boom stiffness has been found to 
occur when the first symmetric and antisymmetric frequencies are equal. 
This tension yields the lightest boom for a particular lowest deployed natural 
frequency. 

It can be seen in Fig. 3 that the first-mode frequency falls off rela- 
tively sharply below the maximum frequency point, but remains relatively 
flat for total blanket tensions up to approximately 5 0% of the boom-buckling 
load. Because of this phenomena, the blanket tension may be increased 
above the maximum frequency point, thus allowing the user to increase the 
antisymmetric first-mode frequency while only slightly decreasing the sym- 
metric first mode. This tension increase has a minimal affect on the over- 
all weight of the roll-up array and allows the user to operate on the flat part 
of the symmetric frequency curve, where variations in the blanket tension 
have little effect on the array frequency. 

The program RUSAP has provisions for calculating the performance 
parameters for any selected boom and blanket tension, and, in addition, has 
provisions for calculating the optimum blanket tension. The optimum blanket 
tension is defined by the condition that symmetric and antisymmetric fre- 
quencies are equal. Each deployed array analysis requires the calculation of 
the natural frequencies of the combined tensioned-blanket/axial-loaded boom 
system that makes up the deployed array. In addition to the optimum tension 
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analysis, the program can also be used to determine the boom stiffness that 
results in the lowest deployed frequency being equal to a desired minimum 
deployed natural frequency. 

The modal analysis of the deployed solar array is based on a multiple- 
degr ee-of-fr eedom, finite-element representation of the boom, beam, and 
blanket components, as shown in Fig. 4. Each blanket is modeled by 10 
rectangular finite elements, which describe the out-of-plane stiffness caused 
by the imposed blanket tension. The boom and beam components are simi- 
larly modeled by standard beam-column elements, which describe both the 
bending stiffness and the geometric stiffness caused by the axial preload. 
Although the blanket bending stiffness is neglected in the above analysis, 
more sophisticated 250-degree-of-freedom analyses, which include the bend- 
ing stiffness, have shown that the simple model leads to less than 1% errors 
in the first-mode frequencies. 

The stiffness and mass matrices for the overall array are developed 
by combining the element stiffness and consistent mass matrices for the ele- 
ments described above. The generation procedure allows for the following 
arbitrary parameters: 

(1 ) Array length. 

(2) Array width. 

(3) Blanket weight/unit area. 

(4) Boom weight/unit length. 

(5) Beam weight/unit length. 

(6) Boom stiffness. 

(7) Beam stiffness. 

(8) Blanket tension. 

After initialization of the above parameters by the root-finding routines, 

the natural frequencies of the array are determined by solution of the usual 
eigenvalue problem with a very fast eigenvalue algorithm. 
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B. Support Structure Analysis 


The support structure analysis uses scaling equations to extrapolate 

the size and weight of the support- structure components of a proposed array 

from those of a prototype array. The configuration used as the baseline for 

2 2 

this analysis is the 23 m (250 ft ) engineering prototype described in Ref. 1. 
Essentially, a dimensional-analysis approach was taken to determine the 
scale factors that are applied to the structural elements when design condi- 
tions change. The total weight of the reference array is broken down into 
12 components: 


W. = w u i+ w ,+W , + W, , + W ,+W ,+W , + W L ,+W ,+W , + W,, +W 

1 bol si cl bl ssl esl csl bel al nsl bktl negl 


( 1 ) 


where 


W 


= weight of boom (determined by modal analysis program) 


W 


= weight of storage drum shell 


W j = weight of end caps on storage drum shell 


W ^ = weight of bearings 


W 


s s 


1 = weight of support shaft 


W 


es 


1 = weight of end supports 


= weight of center support 


W 


= weight of beam (leading-edge member) 


W 


j = weight of boom actuator 


W 


= weight of slip rings and harness 


Wfakti = weight of solar-array blanket 

W 1 = weight of NEG'ATOR spring-mechanism hardware 
° (constant-force spring) 


(Subscript 1 refers to the prototype reference array; the alphabetic 
subscripts refer to the structural components. ) 
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The weight of a second array — differing from the reference array in 
geometry, material properties, and applied inertial loads — can be broken 
down in a similar manner: 


W, = W. , + W ,+'W , + W,,+ W , + W , + W ,+ W, , + W , + W , + W ul +W 
2 bo2 s2 c2 b2 ss2 es2 cs2 be2 a2 ns2 bkt2 neg2 


( 2 ) 


The ratio = W^/W^ can be written as 


W, W u ,W, , W , W , W , W . W, , W, . W , W . 

w _ 2 _ bo2 bol , s2 s 1 c 2 c 1 b2 bl ss2 ssl 

t “ W. _ W. ,W, W . W, W , W. W. . W. W , W. 

1 bol 1 si 1 cl 1 bl 1 ssl 1 


ww ww w w ww ww 

, es2 esl cs2 csl be2 bel a2 al ns2 nsl 

T “JIT -TIT TIT TIT • TTT itr • TIT tit * 


w 


esl 


w. 


1 


w 


csl 


w. 


W, 


bel 


W, 


1 


w , w. 

al 1 


W 


ns 1 


W. 


1 


w w w w 

, bkt2 bktl , neg2 negl 

W u . W, + W W. 

bktl 1 negl 1 


(3) 


where W, . /W, , W , /W, , W , /W, , etc. , are the fractions of the total weight 
bol 1’ si' 1’ cl' 1’ * 6 

of the reference array contributed by the individual components. A simplified 

notation is used to refer to the ratios of the weights of the components, 

W u 0 /W u = W, , W ,/w , = W , W ,/W = W , etc. 

bo2 / bol bo’ s2 ' si s’ c2' cl c’ 

Geometric scale factors are introduced to define the change in size of 
the components, and relationships are then established between these geo- 
metric scale factors and the ratio of stresses in the components. The scale 
factors used are as follows: 


\ 


so 




X 

v 


structural section overall dimension scale factor 
structural section material thickness scale factor 
blanket width scale factor 
array width scale factor 
array length scale factor 
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X^. = ratio of blanket tension of array 2 to that of array 1 


X = ratio of launch- acceleration loading of array 2 to that of 
acc 

array 1 


X 


dia 


ratio of boom diameter of array 2 to that of array 1 


X = ratio of total diameter of combined storage drum shell and 
ves ° 

rolled blanket of array 2 to that of array 1 


Because only one acceleration load factor X relates the acceleration 
loadings in the three orthogonal directions (parallel to the storage drum, 
perpendicular to the storage drum and in the plane of the blanket, and per- 
pendicular to the storage drum and normal to the plane of the blanket), this 
factor is an average of these three orthogonal acceleration factors. 

By expressing the weight ratios in terms of material densities and 
volumes, and by expressing the volumes in terms of the scale factors, 

Eq. (3) becomes 


W 


w t " 


bol 


W. 


w 


bo 


W 

w. 


si 


p s2 

(\ so A St A J 

W . 
+ • — — 


P sl 

\ s s h/ 

w i 

L p d ' s Ac ; J 


w bl 

P b2 j 

K:f 

W ssl 

Pss2 / 

W 1 

P 

L bl 

+ *1 

L P ssl * 


fx so 

\ ss 


X st 

ss 



+ 


+ 


+ 


W es 1 

p es2 

l\ so A st X \ 

+ W csl 

P cs2 / x Bt . \ 

W 1 

_ p es 1 

1A A A. | 

\ es es ves/ 

w, 

P cg i \ dia cs ves) 


W bel 

P be2 Lso \ st ^ ' \ 

+ W a! 

p a2 

fx x st x 1//z ) 

w i 

p, , 'be oe h) 
L bel 

W J 

_ p al 

l dia a v / 


W nsl 

Pns2 /x’X ) 

W bktl 

P bkt2 x . 

W 1 

p nsl ' h V [ 

+ w i 

L P bktl h V 


W , 

negl 


Pneg2 

_ p negl 



(4) 
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where p is the material density of the components, the subscripts on the 
geometric scale factors and densities indicate the structural components to 
which they refer, and the superscripts st and so indicate whether the scale 
factor refers to thickness or to overall size, respectively. 

Relationships are then e stabli shed between the geometric scale factors 
in Eq. (4) and the ratios of stresses of the components of the structure. 

This was done individually for each of the components sized by stress and 
load levels, and the results substituted into Eq. (4). Some components t 
however, are not sized by stress and load levels. 

The ratio of the boom weights is determined by the modal analysis 

program. The NEG'ATOR spring-mechanism weight is a function of blanket 
tension. Blanket weight is a function of X^ and X . Slip ring and harness 
hardware is also assumed to be a function of X^ and X , and bearing size is 
a function of support- shaft size. 

The derivation of the relationships for the storage drum shell will be 
considered as a typical example. For inertial loading, beam-bending 
moments are related by 


-r— = X X.W. 

acc h bs 


(5) 


where + ^ s 2^^bktl + ^sl^’ corres P on di n g ra tio of 

shell maximum bending stresses is 


cr 


s2 


s 1 


M 2 C 2 I 1 

M 1 C 1 I 2 


( 6 ) 


where the shell- section moments of inertia are related by 



The ratio of the maximum bending stresses becomes 

cr 0 X X.W, 
s2 _ acc h bs 

%1 K°f M 


(7) 


( 8 ) 
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Critical buckling stresses for the shell are related by 


^2 

^sbl 



( 9 ) 


Two additional relationships are assumed by the requirement that 
buckling and bending stress ratios be equally critical and that the bending- 
stress ratio be a function of the yield- strength ratio, 


and 


o’ , _ c r _ 

sb2 s2 

cr 0- 

sbl si 


cr 

s2 



K 


cr 


s cr 


ys2 
ys 1 


(10) 


( 11 ) 


where a-y^ I s the yield strength of the shell for array 2, o-y g j is that for 

array 1, and K is the ratio of the factors of safety of the two arrays. These 
s . 

St so 

equations (5 through 11) are then solved for the unknowns k , X. , and W, 

s s 

and substituted into Eq. (4). 


bs’ 


Table 1 lists the relationships established between the geometric scale 
factors in Eq. (4) and the ratios of the stresses for all of the components, 
where 


W 


bscss 


W 


bkt2 


W 


bktl 


+ w _ + w 0 +w 

s2 c2 ss 


+ w . + w , + w 

si cl s 


W 


ba 


W 


bo2 


W 


bol 


+ W 

Tw 


a2 

al 


W 


W, u , + W_+W_+W _+W 0 +W,_ + W, _ + W _ + W ^ 

bkt2 s2 c2 s s2 a2 bo2 b2 neg2 ns2 

bscssa W,... + W.+W.+W . + W . + W, , + W. , +W . +W , 

bktl si cl ssl al bol bl negl nsl 
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and 


r. 

1 


radius of shell of reference array 


r Q = radius of combined shell and rolled blanket for reference array 


The equations are then solved for the unknowns X so , X s *", X s *", X so , X s *" , X SO 


t ^st \ so 


\ OKJ x St X St 

X be’ be’ X a ’ 


s s c ss ss es 

. , . . _ . r , , and W, . The ratios of the weights 

es cs be be a bs bscss ba 6 

of the individual components are then determined by direct substitution of 
these values into Eq. (4). Table 2 lists the solutions for the above scale 
factors . 


To determine W, , W, and W, , three additional equations (12 

through 14) were written and solved for these three unknowns. To solve 
these equations, the scale- factor solutions must first be substituted for the 
scale factors so that the only unknown in each equation is either 


W, , or W, : 
bscss ba 


W 


bs W 


w bkti / 

w bsi V 


X h X v p 


bkt2\ 


bktly 


\ W ,1 

P s2 A SO x St x \ 

) w bsl 

P s i ( s s h ) 


( 12 ) 


where is the weight of the blanket and storage drum shell of the refer- 

ence array; 


W 


W 


bs 1 


W 


bscss W 


bscss 1 
W ssl 


(W U J + 


cl 


bs' W 


bscs s 1 




W 


bscss 1 


ss2 L so' 


s s 1 


( 0(0 


(13) 


where W, , is the weight of the blanket, storage drum shell, end caps, 

DS C S S X 

and support shaft of the reference array; and 


W u i W . /P 0 

, ir _ bol ,„ T , al / a2 , x st\ 

W, ~ TTr W, ) + r~i I \ X X 

ba W, , bo W u . , \p , dia a v 


bal 


bal v al 


l ' 2 ) 


(14) 


where W^ a j is the weight of the boom and the actuator of the reference array. 
Once these equations have been solved for the W's, and all scale factors 
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have been determined, the final weight is calculated and the power -to -weight 
ratio is obtained. 

This technique for analyzing the size and weight of a structure has 
inherent limitations. As in all parametric studies, arbitrary decisions have 
to be made as to how to describe the relationships between load, stress, and 
member size. Because of these limitations, the results obtained from the 
support-structure analysis are not intended as a substitute for a complete and 
detailed structural analysis, but as a good initial estimate of the sizing and 
weight of an array (given a required array area). 

III. PROGRAM USAGE 
A. Selection of Program Levels 

RUSAP is designed to calculate the design parameters of roll-up solar 
arrays of the single -boom, two-drum, tensioned blanket configuration. The 
calculated parameters are power-to-weight ratio, blanket tension, boom 
stiffness, structural member weights, first-mode resonant frequencies, 
actuator width, and required array length and width. 

The computer program has three levels of operation, which are listed 
below. Block diagrams for the three levels are shown in Figs. 5, 6, and 7. 

Level 1: After the user inputs the desired first-mode frequency and 

all other required input parameters, the program calculates 
the optimum boom stiffness and optimum blanket tension that 
gives the required first-mode frequency. It then calculates 
weights and sizes of all base structural members. 

Level 2: After the user inputs the desired boom stiffness and all other 

required input parameters, the program calculates the opti- 
mum blanket tension that maximizes the first-mode frequency 
of the deployed array. It then calculates weights and sizes 
of all base structural members. 

Level 3: After the user inputs the desired boom stiffness, blanket 

tension, and all other required input parameters, the pro- 
gram calculates the deployed symmetric and antisymmetric 
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first-mode frequencies. It then calculates weights and sizes 
of all base structural members. 

There is a major difference between Level 1 and Levels 2 and 3. When 
operating Level 1, the user specifies the required frequency, and the program 
picks many boom stiffnesses in the process of establishing the optimum stiff- 
ness for the required natural frequency. Each time the program tries a new 
boom, it must also determine the optimum tension for that boom. Therefore, 
each boom stiffness iteration during Level 1 is equivalent to running one com- 
plete case at Level 2. If the user does not need an exact frequency, he can 
save computer time by running Level 2 after estimating boom stiffnesses that 
will bound his desired frequency. The user can determine these boom stiff- 
ness estimates by referring to Figs. A-9 through A-12 of Ref. 2. However, 
if the user does desire an exact natural frequency, he can save computer 
time when using Level 1 if he makes as accurate an estimate of the optimum 
boom stiffness as possible using the data in Ref. 2. 

Using the UNIVAC 1108 Exec 8 computer system, RUSAP requires 
approximately 20 seconds to execute Level 1, 5 seconds to execute Level 2, 
and 1 second to execute Level 3. 

B. Minimum Tension and Boom Buckling Considerations 

When operating Level 1, the program establishes the smallest boom 
compatible with the natural frequency requirement when the first-mode fre- 
quency is maximized by using the optimum blanket tension. When the 
required natural frequency is very low, the optimum blanket tension may be 
lower than the minimum tension required to roll the blanket on the drum dur- 
ing retraction. The minimum tension for conventional substrates with bonded- 
on solar cells is defined approximately by the following empirical relationship: 


tension per blanket 
blanket width ~ 


7.2 N/m (0.5 lb /ft) 


A warning message is printed in the program output whenever the above 
criterion is not met. 

If a minimum -weight array is desired and the optimum tension is below 
the minimum tension, the user has a number of options. One is to develop a 
more flexible solar-cell blanket or to add a positive motor drive for array 
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retraction. Another option is to set the blanket tension equal to the minimum 
and then to increase the size of the boom until the boom is not in danger of 
buckling. The critical buckling load for the boom is derived in Appendix C 
of Ref. 3 and is given by: 


tt 2 EI 


critical 


where the symbols and units are the usual ones for column buckling. Recall 
that a plot of first-mode frequency vs percent of buckling load is given 
earlier in Fig. 3. When a particular tension is selected, the user can obtain 
the array that operates at the desired percent of buckling by using the 
LEVEL = -3 option of the program with the PCTBKL parameter negative. 

In this case the program uses the following expression to define the boom El: 


boom El = 


200 X (blanket tension) X (array length) 
2 

tt X (percent buckling) 


2 


The LEVEL = -3 option is also useful when the user wishes to operate 
at a tension somewhat higher than the optimum tension to reduce the sensi- 
tivity of the first-mode frequency to small changes in blanket tension. This 
philosophy is discussed earlier in Subsection II-A. After running Level 1 to 
determine the optimum boom, the user can use LEVEL = -3 with the 
PCTBKL parameter positive to obtain the array with the tension given by: 


blanket tension 


2 

TT 


X (boom El) X (percent buckling) 

^ 

200 X (array length) 


Note that the program always prints the boom buckling load and the percent 
of buckling operating point in the output listing of deployed array parameters. 

C. Array Geometry Input Options 

In all levels of the program there are two ways of inputing the geometric 
requirements of the solar array. The first allows the user to specify the 
total blanket area and array width, whereas the second allows the user to 
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specify the width and length of the individual blankets. The two requirement 
formats are primarily associated with the stage at which the program is 
being applied in the design process. 

The total blanket area and array width are suited for use early in the 
design process when the primary requirements are configuration constraints 
on the array width and desired power output, which implies a total blanket 
area. Later in the design process specific solar-cell module and bus -bar 
layouts lead to detailed specifications on the dimensions of the individual 
blankets. At this point the second format is more suitable. 

With both of the geometric requirement input formats it is necessary 
for the program to determine some of the array dimensions. This is because 
the size of the boom affects the allowable spacing between the blankets. When 
the blanket area and array width are specified, the boom size determines the 
blanket width and the array length, whereas when the blanket length and 
blanket width are specified, the boom size determines the array width. 

Since the actuator width, and therefore the width between blankets, is 
a function of the boom diameter, the program calculates an appropriate 
actuator width for each boom stiffness. The equation used is: 


actuator width 


/ proposed boom diameter \ / prototype 

\ prototype boom diameter / \ actuator width 


When the user inputs the total blanket area and the array width, the 
program calculates the blanket width and array length using 


blanket width 


(array width - actuator width) 
2 


array length = [ a , re ^ :-j7r] + blanket leader 

7 & [array width - actuator widthj 

In the above expression the blanket leader is the total additional length 
of the substrate material extending beyond the celled blanket, and used to 
connect the blanket to the drum and leading edge member (Fig. 1). 
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When the user inputs the blanket width and blanket length, the equations 
used for calculating array width and array length are: 

array width = 2 X (blanket width) + actuator width 
array length = blanket length + blanket leader 

When choosing the length and width for a proposed array it should be 
noted that the ratio of length to width, or aspect ratio of the array, influences 
the array weight. In an earlier study (Ref. 2), the optimum aspect ratio was 
found to be primarily determined by the required lowest natural frequency of 
the array, and essentially independent of the array size. 

In addition to the aspect ratio and first mode natural frequency there 
are many other parameters that affect the weight and cost of a proposed 
array. Many of these parameters may be varied by the user and are described 
in the next section. 

IV. DATA PREPARATION 

RUSAP is a main program designed to be used in batch mode. The 
input parameters that define the roll-up array's size, material properties, 
and desired first-mode natural frequency are read in from data cards using 
namelist names DATA1 and DATA2. Any number of array configurations 
can be generated in a single computer run by including one or more data sets 
in the run stream. A da.ta set consists of a DATA1 group and its associated 
DATA2 group(s). A data group is the collection of cards that define the 
parameters associated with namelist name. 

The FORTRAN name rule, by which any symbolic name that begins 
with I, J, K, L,, M, or N is automatically an integer, is assumed. 

A. NAMELiIST/DATAI General Parameters 

NAMELIST/DATAl/ contains twenty parameters that may be assigned 
values in any DATA1 group.- All of the parameters must be assigned a value 
in the first DATA1 group. Since the parameter assignments of previous data 
groups are used if not reassigned, DATA1 groups following the first need 
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contain only those parameter assignments that are changes from the previous 
parameter assignments. However, it should be noted that the program auto- 
matically updates the parameters BKTL, BKTW, AREA, WIDTH, EIBOM1, 
and TENS at the end of each case by setting them equal to the "Deployed 
Array Parameters," which are printed in the output. These updated 
parameter values are used in the following case unless the user specifies 
new values. j 

The following parameters belong to NAMELIST/DA TAl/ and may be 
included in any order. Data may be in either SI or English units as selected 
using the parameter NUNITS. 

NSETS: an integer to be used in the first DATA1 group to indicate the 

number of cases to be run. RUSAP will read NSETS data 
sets including the first. 

NUNITS: an integer used to define the input units: 0 for English units 

and 1 for SI units. 

LEVEL: an integer used to set the level of operation of the program. 

There are three levels available to the user (1, 2, 3). The 
operation of the three levels is described in Subsection III-A. 
For Level 3 the sign of LEVEL is used to select between 
input parameter options. For LEVEL = 3 the parameters 
EIBOM1 and TENS are used to define the array. For 
LEVEL = -3 the array is defined by the parameters. PCTBKL 
and either EIBOM1 or TENS depending on the sign of PCTBKL. 

POWER: the required power output of the roll-up array; used to calcu- 

late the specific power in W/kg (W/lb). (Units are watts. ) 

FREQR: the required first-mode natural frequency. (Units are Hertz.) 

(Level 1 only. ) 

NTYPE: an integer used to select the array fixed dimensions. If 

NTYPE = 0, BKTL and BKTW are used. If NTYPE = 1, 
AREA and WIDTH are used. 

BKTL: the required blanket length. (Units are meters or feet. ) 

(Required when NTYPE = 0. ) 

BKTW: the required blanket width for each blanket. (Units are 

meters or feet. ) (Used when NTYPE = 0. ) 
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AREA: 


WIDTH: 
TLEADR: 
PCTBKL: 

TENS: 
EIBOM1 : 

EBOOM: 

RBOOM: 

EIEFF: 


2 

the total blanket area of the array. (Units are meters or 
feet^. ) (Used when NTYPE = 1. ) 

the total width of the array. (Units are meters or feet. ) 
(Used when NTYPE = 1. ) 

the total inboard plus outboard blanket leader length. Typi- 
cal value is 0.3 m (1.0 ft). (Units are meters or feet. ) 

a Level 3 parameter used to obtain an array that operates at 
a given percent of buckling. The absolute value of PCTBKL 
is the desired percent of buckling. The sign of PCTBKL is 
used to select between TENS and EIBOM1 as the other 
deployed array stiffness parameter. When PCTBKL is 
positive the array is defined by PCTBKL and EIBOM1; when 
negative, PCTBKL and TENS are used. (LEVEL = -3 only.) 

the tension in each blanket. (Units are newtons or pounds.) 
(LEVEL = 3 or -3 only. ) 

the stiffness (El) of the deployable boom. For Level 1, this 

input boom stiffness is used as a starting value for the boom 

optimization. For Levels 2 and 3, the stiffness is used as 

2 

the required boom stiffness. (Units are Newton -meters or 
2 

pound-feet . ) 

the deployable boom material modulus of elasticity (E). The 

value for steel is 0.2 X 10 12 N/m 2 (0.4175 X 10 10 lb/ft 2 ). 

2 2 

(Units are Newtons /meter or pounds/foot . ) 

the deployable boom material density. The value for steel 
is 7930 kg/m 3 (495 lb/ft 3 ). (Units are kilograms /meter 3 or 
pounds / foot . ) 

the structural efficiency of the deployable boom. Since the 
BI-STEM boom is split into two wrapped, open-tube sec- 
tions, the effective section moment of inertia (I) is less than 
that of a closed tube with the same cross-section dimensions. 
Typical value is 0.80 efficiency. (Unitless quantity. ) 
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DPT: 


the diameter of the deployable boom divided by the boom 
material thickness. The boom wall thickness is twice the 
material thickness. Typical value is 200. (Unitless 
quantity. ) 

ROARAY; the average weight per unit area of the array blankets. 

Typical value is 0. 878 kg/m^ (0. 18 lb/ft^). (Unit s are 
Z Z 

kilograms/meter or pounds/foot . ) 

NDATA2: an integer that specifies how many DATA2 groups are 

included in the data set. 0 causes no base structure to be 
analyzed, 1 reads one DATA2 group, 2 reads two DATA2 
groups, etc. 

B. NAMELIST/DATA2 General Parameters 

NAMELIST/DA TA2/ contains the following 19 parameters that may be 
included in any DATA2 group in any order. They are all ratios of the array 
component properties in the proposed array to those in the prototype array. 
These parameters are all assigned nominal values of 1 . 00 prior to reading 
the first DATA2 group in the first data set. Since the parameter assignments 
of previous DATA2 groups are used if not reassigned, DATA2 groups follow- 
ing the first need contain only those parameter assignments that are changes 
from those in the previous DATA2 group. This is true even when the pre- 
vious DATA2 group is in a previous data set. Note, however, that the 
number of DATA2 groups in each data set must equal NDATA2, even if the 
DATA2 parameters are a repeat of those in the previous data set. 

the ratio of the material density of the storage drum in the 
proposed array to that in the prototype array. 

the ratio of the modulus of elasticity of the storage drum in 
the proposed array to that in the prototype array. 

the ratio of the yield strength stress of the storage drum in 
the proposed array to that in the prototype array. 

the ratio of the material density of the storage drum end 
caps in the proposed array to that in the prototype array. 

the ratio of the modulus of elasticity of the storage drum 
end caps in the proposed array to that in the prototype array. 
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RHOS: 

ES: 

SIGMS: 

RHOC: 

EC; 



SIGMC: 

RHOSS: 

ESS:. 

SIGMSS 

RHOBE: 

EBE: 

SIGMBE: 

RHOES: 

EES: 

SIGMES: 

RHOCS: 

ECS: 

SIGMCS : 

VLAMAC 


the ratio of the yield strength stress of the storage drum end 
caps in the proposed array to that in the-prototype array. 

the ratio of the material density of the drum support shaft in 
the proposed array to -that in the prototype array. 

the ratio of the modulus of elasticity of the drum support 
shaft in the proposed array to that in the prototype array. 

the ratio of the yield strength stress of the drum support 
shaft in the proposed array to that in the prototype array. 

the ratio of the material density of the leading edge beam in 
the proposed array to that in the prototype array. 

the ratio of the modulus of elasticity of the leading edge beam 
in the proposed array to that in the prototype array. 

the ratio of the yield strength stress of the leading edge beam 
in the proposed array to that in the prototype array. 

the ratio of the material density of the drum end supports in 
the proposed array to that in the prototype array. 

the ratio of the modulus of elasticity of the drum end sup- 
ports in the proposed array to that in the prototype array. 

the ratio of the yield strength stress of the drum end sup- 
ports in the proposed array to that in the prototype array. 

the ratio. of the material density of the center support in the 
proposed array to that in the prototype array. 

the ratio, of the modulus of elasticity of the center support 
in the proposed array to that in the prototype array. 

the ratio of the yield strength stress of the center support 
in the proposed array, to that in the prototype array. 

a factor used to scale the launch vibration level from that 
set within the program. The factor A. , the ratio of the 
launch acceleration loading of the proposed array to the load- 
ing used in the prototype array, is used in calculating array 


18 


JPL Technical Memorandum 33-634 



stresses and member sizes (Subsection II-B). The factor 


\ is defined by: 
acc 1 


\ 

acc 


VLAMAC X 


( Total weight of 
prototype array 
Total weight of 
proposed array 


i/2 


The above square root factor implements a general rule of 
thumb that the launch acceleration loading for a subsystem on 
a given launch vehicle varies with the square root of the sub- 
system's mass. In other words, a small light subsystem 
sees a higher excitation level than a large heavy one. on the 
same spacecraft. VLAMAC is used to change this relation- 
ship to account for different acceleration levels associated 
with different spacecraft/launch-vehicle combinations. The 
levels used on the reference array are described in Figs. 8 
and 9. 

The material properties of the individual components in the prototype 
array are listed in Table 3. They are to be used when calculating the param- 
eter ratios previously described, 


V. EXAMPLES AND SAMPLE OUTPUT 


To demonstrate the program's use, three example problems are con- 
sidered in Table 4. Note that Case 1 has two $DATA2 cards associated with 
two different base structure analyses. 

The following data cards were used to input the above cases: 

$DATA1 NSETS=3,NUNITS=1,LEVEL=1,POWER=2500.,FREQR=.077,NTYPE=1,AREA=23.225, 
WIDTH=2.516,TLEADR=.305,EIBOM1=2000.,EBOOM=0.2E+12,RBOOM=7930.,EIEFF=.80, 
DPT=200 . , ROARAY= . 878, NDATA2=2$ 

$DATA2 VLAMAC=. 5$ 

$DATA2 VI_AMAC= 1 . $ 

$DATA1 LEVEL=2,EIBOM1=1500.,NDATA2=1$ 

$DATA2 VLAMAC=1.$ 

$DATA1 NUNITS=0, LEVEL=-3, POWER=10000 . ,NTYPE=0, BKTL=8 1 . 46, BKTW=6 . 16,TLEADR=1 . 
PCTBKL=-50 . , TENS=3 . 25, EBOOM= . 418E+10, RBOOM=495 . , ROARAY= . 18$ 

$DATA2 VLAMAC=1.$ 
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A. 


Description of Example Cases 


Case 1 is an example of Level 1. The program determines the mini- 
mum weight design for the given size and natural frequency. Case 2 is an 
example of Level 2. The program uses the desired boom stiffness and deter- 
mines the blanket tension that maximizes the frequency. Case 3 is an example 
of Level 3. The user specifies the PCTBKL and blanket tension, and the 
program calculated the boom stiffness and the first mode symmetric and 
antisymmetric frequency. This run is typical of a design with a configuration 
for which the optimum boom requires an unacceptably low blanket tension. In 
other words, the optimum blanket tension is too low to roll up the blanket 
during array retraction. A minimum weight design is obtained by setting the 
tension equal to the minimum recommended tension of 7.2 N/m (0.5 lb/ft) 
and increasing the boom stiffness until the boom is not in danger of buckling 
(Subsection III-B). In this case the boom is selected to operate at 50% of its 
buckling load. 

B. Typical Level 1 Printout 

******* ************************************* *************************** ********* 
******************************************************************************** 

***** POWERS 2500. WATTS AREA- 23.2 M**2 WIDTH: 2.516 M ***** 

***** REQUIRED FREQUENCY: .077 HZ ***** 

******************************************************************************** 
******************************************************************************** 


****INPUT PARAMETERS FOR 0EPL0YE0 ROLL-UP SOLAR ARRAY**** 


REQUIRED POWER 

= 

2500. 

WATTS 




BLANKET AREA 

= 

250.0 

FT**2 

( 

23.22 

M**2> 

ARRAY WIDTH 

= 

8.25 

FT 

( 

2.516 

M> 

leader length 

= 

1.00 

FT 

< 

.305 

M> 

BLANKET UNIT WT. 

= 

• 160 

L8/FT **2 

( 

.88 

KG/M**2! 

REQUIRED FREQ 

= 

.077 

HZ 




Initial boom (ei> 

= 

4837. 

LB FT **Z 

( 

2000. 

NM**2) 

BOOM MODULUS (E) 

= 

.4175+10 

LB/FT **2 

( 

.200+12 

N/M**2) 

BOOM MATL DENSITY: 

495.0 

LB/FT **3 

( 

7930. 

KG/M**3 

BOOM EFFICIENCY 

= 

.80 

effective. 

(I) 



DIA TO THICKNESS 

= 

200. 

IN/IN 





♦•♦♦PRINTOUT OF PROGRAM OPTIMIZATION ITERATIONS**** 
(ENGLISH UNITS-L8-FT.SEC.HZ) 


TENSION 

OPTIMIZATION 

ITERATIONS 




I 

TENSION 

BOOM (El) 

FREQ DIF 

STM FREQ 

ASYM FREQ 

1 

.41857+00 

.48368+04 

#29618-01 

.62013-01 

.32395-01 

2 

.30247+01 

.48368+04 

.20164-01 

.10710+00 

.86940-01 

3 

.30247+01 

.48368+04 

.20164-01 

.10710+00 

.86940-01 

4 

.48027+01 

,48366+04 

-.41058-03 

.10862+00 

.10903+00 

5 

,47672+01 

.48368+04 

-.30647-04 

.10873+00 

.10876+00 


END TENSION OPTIMIZATION 
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BOOM STIFFNESS OPTIMIZATION ITERATIONS 
I TENSION eOOM (El) FREQ DIF FREQ FREQ REQ 

0 .47672+01 .48368+04 .31743-01 .10874+00 .77000-01+ 


TENSION 

OPTIMIZATION 

ITERATIONS 





I 

TENSION 

BOOM (El) 

FREQ OIF 


STM FREQ 

ASYM FREQ 

1 

.44867+00 

.51995+04 

.30371-01 


.64153-01 

.33782-01 

2 

.32131+01 

.51995+04 

.21199-01 


.11057+00 

.89373-01 

3 

.32131+01 

.51995+04 

.21199-01 


.11057+00 

.89373-01 

4 

.51614+01 

.51995+04 

-.76944-03 


.11215+00 

.11292+00 

5 

.50932+01 

.51995+04 

.13062-03 


.11237+00 

.11224+00 





end 

TENSION OPTIMIZATION 

TENSION 

OPTIMIZATION 

ITERATIONS 





I 

TENSION 

BOOM (El) 

FREQ OIF 


SYM FREQ 

ASYM FREQ 

1 

.10953+00 

.12092+04 

.15541-01 


.32114-01 

.16573-01 

2 

.79294+00 

.12092+04 

.11377-01 


.56571-01 

.45194-01 

3 

.79294+00 

.12092+04 

.11377-01 


.56571-01 

.45194-01 

4 

.13099+01 

.12092+04 

-.23155-03 


.57842-01 

.58073-01 





END 

TENSION OPTIMIZATION 


BOOM stiffness optimization iterations 

I TENSION BOOM (El) FREQ DIF FREQ FRFQ REQ 

1 .13099+01 .12092+04 -.19043-01 .57957-01 .77000-01+ 


TENSION OPTIMIZATION ITERATIONS 


I 

TENSION 

BOOM (El ) 

FREQ DIF 

SYM FREQ 

ASYM FREQ 

1 

.22751+00 

.25694+04 

.21986-01 

.46005-01 

.24019-01 

2 

.16370+01 

.25694+04 

.15881-01 

.80327-01 

.64446-01 

3 

.16370+01 

.25694+04 

.15881-01 

.80327-01 

.64446-01 

4 

.26803+01 

.25694+04 

-.26172-03 

.81979-01 

.82241-01 


ENO TENSION OPTIMIZATION 


BOOM STIFFNESS OPTIMIZATION ITERATIONS 

1 TENSION BOOM (El) FREQ DIF FREQ FRrQ REQ 

2 .26803+01 .25694 + 04 .51103-02 .82110-01 .77f100-01« 


TENSION 

optimization 

ITERATIONS 




I 

TENSION 

300M (El) 

FREQ DIF 

SYM FREQ 

ASYM FREQ 

1 

.20282+00 

.22816+04 

.20571-01 

.43481-01 

.22910-01 

2 

.14627+01 

.22816+04 

.15045-01 

.76101-01 

.61056-01 

3 

. 14627 +01 

.22816+04 

.15045-01 ✓-N. 
-.43306-08 

end 

.76101-01 

.61056-01 

4 

•23931+01 

.22816+04 

.77622-01 .77622-01 

TENSION OPTIMIZATION 


BOOM STIFFNESS OPTIMIZATION ITERATIONS 
I TENSION BOOM (El) frREQ DIF FREQ FREQ REQ 

3 .23931 + 01 .22816+04 .52191-03Q .77622-01 (3) .77000-01* 


There are many error messages that can appear in the above section. 
However, they may not affect the final answer. Without explaining the 
details of all the error messages, an easy check can be made to determine 
whether the program has run ^successfully. If the frequency differences 
marked (T) and ( 2 ) are less than approximately 5% of the calculated fre- 
quency (T) then both the sections have run successfully. If both of these 
sections have been completed successfully, the optimization has been success- 
ful, even if there are error messages. If the run is not successful, check 
your input data. If the data is correct, your configuration may have caused 
the eigenvalue routine to fail; and the only thing that can be done is to try a 
different boom stiffness or required frequency and see if that helps 
optimization. 
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****»*****OEPLOYED ARRAY PARAMETERS*** ******* . . 

REQUIRED POWER = 2500. WATTS 

BLANKET AREA ' - 25U.0 FT**2 ( 23.22 M**2> 

BLANKET. WIDTH = 3.67 FT. ( 1.160 M) 

BLANKET LENGTH = 32.32 FT ( 9.85B M) 

ARRAY WIDTH = 6.25 FT < 2.516 Ml 

ARRAY LENGTH = 33.32 FT ( 10.163 M) 

ACTUATOR WIDTH = '•' .51 FT ( .157 M) 

LEADER LENGTH = 1.00 FT ( .305 M) 

BLANKET UNIT WT. = .180 LB/FT**2 ( .88 KG/M**2) 

BLANKET TENSION = 2.39 LB/BLANKET l 10.6 N/BLANKET) 

APPLIED BOOM . LOAD= .23 OF BOOM BUCKLING LOAD 

BUCKLING LOAD = 20.26 LB ( 90.2 N) 

BOOM STIFFNESS EI= 2262. LB FT**2 ( 943. NM**2> 

BOOM DIAMETER . = . : 1.38 IN- < .0350 M) 

BOOM THICKNESS = .007 IN .( . 00018. M) 

BOOM UNIT WEIGHT = .19 LB/FT < .29 KG/M) 

BOOM MODULUS (El = .4175+10 LB/FT**2 < .200+12 N/M**2> 

BOOM MATL DENSITY^ 495.0 LB/FT**3 < 7930. KG/M**3> 

BOOM EFFICIENCY : .80 EFFECTIVE (I) 

FREQUENCY = .078 HZ. 

• i ; *-. • 

****6ASE STRUCTURE INPUT PARAMETERS**** 


RHOS 

ES 

sigms 

RHOC 

EC 

SIGMC 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

RHOSS 

ESS 

sigmss 

RHObE 

ebe 

SlGMBE 

i.opo 

1.000 

1.000 

l.obo 

1.000 

- l.ooo 

RHOES 

EES 

SIGMES 

RHOCS 

ECS 

SIGMCS 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 


VLAMAC 

.500 



*******CALCULATE0 BASE STRUCTURE PARAMETERS******* 

Q (LAUNCH ACCELERATION^ .50*BASELINE LEVEL) 

TOTAL BLANKET WEIGHT = 46.37 LBS < 21.05 KG) 

DEPLOYED BOOM WEIGHT - 6.45 LBS < 2.93 KG) 

LEADING EDGE BEAM WEIGHT = .74 LBS < .34 KG) 

OUTBOARD END SUPPORT ASSY WT*= 1-89 LBS < .86 KG) 

DRUM SHELL WT. = 3.51 LBS ( 1.59 KG) 

END CAP+GUIOE WT. = 1.68 LBS ( .76 KG) - 

TOTAL DRUM ASSY WEIGHT = 5.19 LBS < 2.36 KG) 

CENTER SUPPORT WT. = .72 LBS ( .33 KG) 

BOOM ACTUATOR WT. = 3.47 LBS ( 1.58 KG) 

BEARING ASSY WT.- = . 1.47 LBS . < .67 KG) , . . 

NEGATOR ASSY WT. = 1.68 LBS ( .76 KG) 

SLIPRING+HARNESS WT= 2.35 LBS ( 1.07 KG) 

SUPPORT SHAFT WT. = .71 LBS ( .32 KG) 

TOTAL CENTER SUPPORT ASSY WT = 10.40 LBS < 4.72 KG) 

TOTAL SOLAR ARRAY WEIGHT = 71.06 LBS < 32.26 KG) 

POWER TO WEIGHT EFFICIENCY = 35.18 WATTS/LB ( 77.49 WATTS/KG) 
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If there is an error message at the point marked (4) in the program, 
it means one or more of the base structure equation solutions did not con- 
verge. If this happens, check all input data for errors. If there are no 
errors, the configuration may be too large or too small for the analysis to 
converge. 


****BASE STRUCTURE INPUT PARAMETERS**** 


RHOS 

ES 

SIGMS 

1.000 

1.000 

1.000 

RHOSS 

ESS 

SIGMSS 

1.000 

1.000 

1.000 

RHOES 

EES 

SIGMES 

1.000 

1.000 

1.000 

VLAMAC 

1.000 




RHOC 

EC 

SIGMC 

1.000 

1.000 

1.000 

RHOBE 

EBE 

sigmbe 

1.000 

1.000 

1.000 

RHOCS 

ECS 

SIGMCS 

1.000 

1.000 

1.000 


♦♦♦♦♦♦♦CALCULATED BASE STRUCTURE PARAMETERS******* (T) 
(LAUNCH ACCELERATION^ 1 . 00 *B ASELINE LEVEL) ^ 


TOTAL BLANKET WEIGHT = 46.37 LBS < 
DEPLOYED BOOM WEIGHT = 6.45 LBS ( 
LEADING EDGE BEAM WEIGHT = .93 LBS ( 
OUTBOARD END SUPPORT ASSY WT.= 4.08 LBS ( 


DRUM SHELL wT. = 

END CAP+GUIDE WT. = 

TOTAL DRUM ASSY WEIGHT 

CENTER SUPPORT WT. = 

BOOM ACTUATOR WT. = 

BEARING ASSY WT. = 

NEGATOR ASSY WT. = 

SLIPRING+HARNESS WT= 

SUPPORT SHAFT WT, = 

TOTAL CENTER SUPPORT ASSY WT = 

TOTAL SOLAR ARRAY WEIGHT = 

POWER TO WEIGHT EFFICIENCY = 


5.74 LBS ( 2.61 KG) 

4.49 LBS ( 2.04 KG) 

= 10.23 LBS < 


1.35 LBS ( 

5.35 LBS ( 
2.49 LBS ( 
1.68 LBS ( 

2.35 LBS ( 
1.21 LBS ( 


14.42 LBS ( 
82.49 LBS ( 


21.05 KG> 
.2.93 KG). 
.. * 42 KG) 
1.85 KG) 

4.65 KG) 


.61 KG) 
2.43 KG) 
1.13 /KG) 
.76 KG) 
1.07 KG) 
.55 KG) 


6.55 KG) 
37.45 KG) 


30.31 WATTS/LB ( 66.76 WATTS/KG) 


TIME FOR RUSAP = 17.07 SEC. 


This section (marked \5j above) is the second base structure analysis 
The deployed array parameters remain the same and new base parameters 
are read. 
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c 


Typical Level 2 Printout 


**************************** ******* ************************ ********************* 
******************************************************************************** 
***** POWER= 2500. WATTS AREA= 23.2 M**2 WIDTH= 2.516 M ***** 

***** 30 OM STIFFNESS= 1500. N*M**2 ***** 

******************************************************************************** 
******************************************************************************** 


****INPUT PARAMETERS FOR DEPLOYED ROLL-UP SOLAR ARRAY**** 


required power 

= 

2500. 

WATTS 



BLANKET AREA 

= 

250.0 

Vt**2 

( 23.22 

M**2) 

array width 

= 

B. 25 

FT 

( 2.516 

M) 

LEADER LENGTH 

= 

1.00 

FT 

( .305 

M) 

BLANKET UNIT WT. 

= 

.180 

LB/FT**2 

( .88 

KG/M**2> 

BOOM STIFFNESS El 

:= 

3628. 

LB FT **2 

( 1500. 

NM**2) 

BOOM MODULUS (E) 

= 

.4175+10 

LB/FT**2 

( .200+12 

N/M**2) 

BOOM MATL DENSITY 

- 

495.0 

LB/FT**3 

( 7930. 

KG/M**3) 

BOOM EFFICIENCY 

= 

.80 

EFFECTIVE (I) ' 


DIA TO THICKNESS 


200. 

IN/IN 




****pr intout of program optimization iterations**** 

(ENGLISH UNITS-LB .FT r SEC »HZ) 


TENSION OPTIMIZATION ITERATIONS 


I 

TENSION 

BOOM (El) 

FREQ DIF 

SYM FREQ 

ASYM FREQ 

1 

.31737+00 

.36276+04 

.25745-01 

.54161-01 

.28416-01 

2 

.22758+01 

.36276+04 

.18272-01 

.94009-01 

.75737-01 

3 

.22758+01 

.36276+04 

.18272-01 

.94009-01 

.75737-01 

4 

.36868+01 

.36276+04 

-.17921-03 

.95726—01 

.95905-01 


END TENSION OPTIMIZATION 


**********deployed array parameters********** 


REQUIRED POWER 

= 

2500. 

WATTS 




BLANKET AREA 

= 

250.0 

FT**2 

( 

23.22 

M**2> 

blanket WIDTH 

= 

3.84 

FT 

( 

1.170 

M ) 

BLANKET LENGTH 

= 

32.59 

FT 

( 

9.939 

M) 

ARRAY WIDTH 

= 

8.25 

FT 

( 

2.516 

M) 

ARRAY LENGTH 

= 

33.59 

FT 

( 

10.244 

M) 

ACTUATOR WIDTH 

= 

.58 

FT 

( 

.176 

M) 

LEADER LENGTH 

= 

1.00 

FT 

( 

.305 

M) 

BLANKET UNIT WT. 

n 

.180 

LB/FT**2 

( 

.88 

KG/M**2) 
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BLANKET TENSION = 

3.69 

LB/BLANKET 

( 

16.4 

N/BLANKET) 

APPLIED BOOM L0AD= 

.23 

OF BOOM BUCKLING LOAD 


BUCKLING LOAD = 

31.71 

LB 

( 

141.1 

N) 

BOOM STIFFNESS EI= 

3628. 

LB FT**2 

( 

1500. 

NM**2 ) 

BOOM DIAMETER = 

1.55 

IN 

( 

.0393 

M) 

BOOM THICKNESS = 

.008 

IN 

( 

.00020 

M) 

BOOM UNIT WEIGHT = 

.24 

LB/FT 

( 

.36 

KG/M) 

BOOM MODULUS (E) = 

.4175+10 

LB/FT**2 

( 

.200+12 

N/M**2> 

BOOM MATL DENSITY:: 

495.0 

LB/FT**3 

( 

7930. 

K6/M**3> 

BOOM EFFICIENCY = 

.80 

EFFECTIVE 

(I) 



FREQUENCY = 

.096 

HZ 





****BASE structure input parameters**** 


RHOS 

ES 

SIGMS 

RHOC 

EC 

SIGMC 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

RHOSS 

ESS 

SIGMSS 

RHOBE 

EBE 

SIGMBE 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

RHOES 

EES 

SIGMES 

RHOCS 

ECS 

SIGMCS 

1.000 

1.000 

1.000 

1.000 

1.000 

i.ooo 


VLAMAC 

1.000 


*******calculated base structure parameters******* 

(LAUNCH ACCELERATION= 1 . 00*BASELINE LEVEL) 


TOTAL blanket weight = 
DEPLOYED BOOM WEIGHT = 
LEADING EDGE BEAM WEIGHT = 
OUTBOARD END SUPPORT ASSY WT.= 


46.36 LBS ( 
8.20 LBS ( 
1.17 LBS ( 
4.07 LBS ( 


21.05 KG) 
3.72 KG) 
.53 KG) 
1.85 KG) 


DRUM SHELL WT. = 5.66 LBS ( 2.57 KG) 

END CAP+GUIDE WT. = 4.47 LBS ( 2.03 KG) 

TOTAL DRUM ASSY WEIGHT = 10.13 LBS ( 4.60 KG) 


CENTER SUPPORT WT. = 
BOOM ACTUATOR WT. = 
BEARING ASSY WT. = 

NEGATOR ASSY WT. = 

SLIPRING+HARNESS WT= 
SUPPORT SHAFT WT. = 
TOTAL CENTER SUPPORT ASSY 


1.46 LBS ( .66 KG) 
6.80 LBS ( 3.09 KG) 

2.47 LBS ( 1.12 KG) 
2.09 LBS ( .95 KG) 
2.37 LBS ( 1.07 KG) 
1.19 LBS < .54 KG) 
WT = 16.37 LBS ( 


7.43 KG) 


TOTAL SOLAR ARRAY WEIGHT 


86.30 LBS (' 39.18 KG) 


POWER TO WEIGHT EFFICIENCY = 28.97 WATTS/LB ( 63.81 WATTS/KG) 


TIME FOR RUSAP = 3.33 SEC. 
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D. Typical Level 3 Printout 


***********************'********************************************************* 
******************************************************************************** 
***** POWER=10000. WATTS BLANKET L= 81.46 FT BLANKET W= 6.160 FT ***** 

***** BOOM STlFFNESS= 8965. LB*FT**2 TENSION: - 3.25 LB ***** 

******************************************************************************** 
******************************************************************************** 


****input parameters for deployed roll-up solar array**** 


required power = 

10000. 

WATTS 




blanket area = 

1003.6 

FT**2 

( 

93.23 

M**2 ) 

blanket width = 

6.16 

FT 

t 

1.879 

M) 

blanket length = 

81.46 

FT 

( 

24.845 

M) 

LEADER LENGTH = 

1.00 

FT 

( 

.305 

M) 

BLANKET UNIT WT. = 

.180 

LB/FT**2 

{ 

.88 

KG/M**2 ) 

BOOM STIFFNESS El: 

8965. 

LB FT**2 

( 

3707. 

NM**2) 

BOOM MODULUS (E) = 

.4180+10 

LB/FT**2 

( 

.200+12 

N/M**2) 

BOOM MATL DENSITY: 

495.0 

LB/FT **3 

( 

7930. 

KG/M**3 ) 

BOOM EFFICIENCY = 

.80 

EFFECTIVE 

(I) 



DIA TO THICKNESS = 

200. 

IN/ IN 




PERCENT BUCKLING : 

50.000 





BLANKET TENSION : 

. 3.25 

LB/BLANKET 

( 

14.5 

N/BLANKET) 

**********OEPLOYED 

ARRAY PARAMETERS********** 


REQUIRED POWER .= 

10000. 

WATTS 




BLANKET AREA = 

1003.6 

FT**2 

( 

93.23 

M**2 ) V 

BLANKET WIDTH = 

6.16 

FT 

( 

1.879 

M) 

blanket LENGTH = 

81 . 46 

"FT 

( 

24.845 

M) . . 

ARRAY WIDTH = 

13.04 

FT 

( 

3.978 

M) 

ARRAY LENGTH = 

82.46 

FT 

( 

25.150- 

Wc 


ACTUATOR WIOTH = 

.72 

FT 

( 

.221 

"Ml 

LEADER LENGTH = 

1.00 

FT 

( 

.305 

M) 

BLANKET UNIT WT. = 

.180 

LB/FT**2 

( 

.88 

KG/M**2) 

BLANKET TENSION = 

3.25 

LB/BLANKET 

( 

14.5 

N/BLANKET) 

APPLIED BOOM LOAD: 

.50 

OF BOOM BUCKLING LOAD 


BUCKLING LOAD : 

13.00 

LB 

( 

57.8 

N) 

BOOM STIFFNESS El: 

8965. 

LB FT **2 

( 

3707. 

NM**2 ) 
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BOOM DIAMETER = 1.94 IN (' ,0493'M) 

BOOM THICKNESS = .010 IN • < .00025 M) 

BOOM UNIT WEIGHT = .38 LB/FT C -.57 KG/M) 

BOOM MODULUS (E) = .4180 + 10 LB/FT**2- ( .200+12 N/M**2) 

BOOM MATL DENSITY= -495.0 LB/FT++3 < -793^0. KG/M**3> 

BOOM EFFICIENCY .. = .80 EFFECTIVE (I) 

CALCULATED SYM. FREQ.= .018 HZ.' CALCULATED ASYM. FREQ.= .029 HZ. 


****BASE STRUCTURE INPUT PARAMETERS**** 


RHOS 

ES 

SIGMS 

RHOC 

EC 

’ SIGMC 

1 .000 

1.000 

1.000 

1.000 

1.000 

i .ooo ' 

RHOSS 

ESS 

SIGMSS 

RHOBE' 

EBE 

SIGMBE 

1 .000 

1.000 

1.000 

1.000 

1.000 

1.000 , 

RHOES 

EES 

SIGMES 

RHOCS 

ECS 

■SIGMCS . 

1.000 

VLAMAC 

1.000 

1.000 

1.000 

1.000 

-1.000 

. i.000 . 


*******CALCULATEO BASE STRUCTURE PARAMETERS******* 

(LAUNCH ACCELERATION^ .50+BASELINE LEVEL) 

TOTAL BLANKET WEIGHT .= 182.86 LBS .( 83.02 KG) 

DEPLOYED BOOM WEIGHT = 31.64 LBS ( 14.37 KG) 

LEADING EDGE BEAM WEIGHT =' '' 2.42 LBS' ( 1.10 KG) 

OUTBOARD END SUPPORT ASSY WT.= 12.37 LBS ( 5.62 KG) 

DRUM SHELL WT. = 19. 50. LBS ( 8.85 KG) 

ENO CAP+GUIDE WT. = 12.84 LBS ( 5.83 KG) 

TOTAL DRUM ASSY WEIGHT = 32.34 LBS ( 14.68 KG) 

CENTER SUPPORT WT. = 4.03 LBS ( 1.83 KG) 

BOOM ACTUATOR WT.. = 17.09 LBS ( 7.76 KG) 

BEARING ASSY WT. = 5.22 LBS ( 2.37 KG) 

NEGATOR ASSY WT. = 1.96 LBS ( .89 KG) 

SLIPRING+HARNESS WT= 9.19 LBS.< 4.17 KG) 

SUPPORT SHAFT WT. = 4.03 LBS ( 1.83 KG) 

TOTAL CENTER SUPPORT ASSY WT = 41.53 LBS < 18.85 KG) 

TOTAL SOLAR ARRAY WEIGHT =; 303.16 LBS ( 137.63 KG) 

POWER TO WEIGHT EFFICIENCY = 32.99 WATTS/LB ( 72.66 WATTS/KG) 


TIME FOR RUSAP = .86 SEC. 
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Table 1. Geometric scale factor and stress-level relationships 






















































































Table 4. Data defining three example cases 


Parameter 


Nominal values preset 
in program 


Case 1 


Case 2 


Case 3 


$DATA1 


NSETS 


3 

_b 


NUNITS 


1 

— 

0 

LEVEL 


1 

2 

-3 

POWER 


2500. 

— 

10000. 

FREQR 


0.077 

* 

* 

NTYPE 


1 

• — ■ 

0 

BKTL 


- *a 

* 

81.46 

BKTW 



* . 

6.16 

AREA 


23.225 

■ — 

❖ 

WIDTH 

NO DATA 1 VALUES 

2.516 

\ — ’ 

* 

TLEADR 

ARE PRESET 

0.305 

— 

1.0 

PCTBKL 


❖ 


-50. 

TENS 


* 

-* 

3.25 

EIBOM1 


2000. 

1 500. ; 

* 

EBOOM 


0.20E+1 2 

• . . - — . , 

0.418E+10 

RBOOM 


7930. 

— 

495. 

EIEFF 


0.80 


— 

DPT 


200. 

— 

— 

ROARAY 


0.878 

— 

0.18 

NDATA2 


2 

.1 



$DATA2 


RHOS 

1.0 





ES 

1.0 

— 

— 

— 

— 

SIGMS 

1.0 

— 

— 

— 

— 

RHOC 

1.0 

— 

— . 

— 

— 

EC 

1.0 

— 

— _ 

— 

— 

SIGMC 

1.0 

— 

— 

— 

— 

RHOSS 

1.0 

— 

— 

— 

— 

ESS 

1.0 

— 

— 

— ' ■ 

— 

SIGMSS 

1.0 

— 

— 

. — 

— 

RHOBE 

1.0 

— 

— 


— 

EBE 

1.0 

— 

— 

• — ' 

— 

SIGMBE 

1.0 

— 

— : 

— 

— 

RHOES 

1.0 


— 

— . 

— 

EES 

1.0 

' — 



— 

SIGMES 

1.0 

— 

— 

— • 

— 

RHOCS 

1.0 

— 

■— 

— 

— 

ECS 

1.0 

— 

— 

— 

— 

SIGMCS 

1.0 

- 

— 

— 

— 

VLAMAC 

1.0 

0.5 

1.0 

; — • 

— 


a (*) indicates that no value need be assigned, since this parameter is not used for this case. 
k(— ) indicates the previous or preset values of the variables are to be used. 
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Fig. 1. Deployed array 


(o) TORSION (ANTISYMMETRIC) 



(b) BENDING (SYMMETRIC) 



Fig. 2. Deployed array mode shapes 
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TENSION, % OF BOOM BUCKLING LOAD 


Fig. 3. Typical plot of frequency 
vs blanket tension 



Fig. 4. Typical finite - element models: 

(a) symmetric; (b) antisymmetric 
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Input Deployed 
Array Parameters 


Calculate Actuator Width 


If NTYPE = 0 Calculate Array 
Width and Array Length 


If NTYPE = 1 Calculate Blanke 
Width and Array Length 


Choose Blanket Tension 


Modal Analysis Program 


Fundamental 
Symmetric Frequency 
Asymmetric Frequenc 


Output Deployed 
Array Analysis Data 


Input Base Structure 
Parameters 


Calculate 

Acceleration Factor 


Calculate Size and 
Weight of Base Structure 


- — ' — (Calculated Weight 
Previous Calculated Weight) 
Less Than 

^SU-lovable Difference^ 


Output Base Structure 
Analysis Data 


u*. 

Yes 

An Additional Base\. 
Analysis to be Run^*" 
^\For This Cas£^^ 

No 


Stop 

Fig. 6. 

Level 2 flow chart 
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10 100 1 70 700 1000 10,000 

FREQUENCY, Hz 


Fig. 8. Qualification sine vibration levels 
for the GE rollup array 



10 100 1000 10,000 


FREQUENCY, Hz 


Fig. 9. Qualification random vibration levels 
for the GE rollup array 
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C RUSAP (MAIN DRIVING PROGRAM) 

C 

C 

REAL L 

COMMON/JSP1/EIEFF .FREQ. DIABOM»THIBOM.ROBOOM.WIDCE»DEGOP»WPANEL» MAX 
HTR.NDATA2.ROARAY»NW»NL.L.WIDTH»FREQR»RBOOM,EBOOM»RPT, POWER. AREA »T 
2LEADR*BKTW»BKTL*NTYPE 

COMMON/ JSP2/GEDI AB. GEW IOC. GEW » GET . GEE I BE » GEROBE 

COMMON/ JSP6/VLAMAC . VLAMT . VLAMH * RHOBE . AK.BE . SIGMBE »EBE » VLMDI A 

COMMON/ J5P7/PCTBKL 

1 FORMAT (5F 10*5. 415) 

2 FORMAT (4110) 

4 FORMAT (1H0. 'BOOM STIFFNESS El='.F10.0»* LB FT**2 C.F8.0.' NM** 

12 ) •) 

5 FORMAT <//lH0» 'FREQUENCY ='.F10.3>' HZ') 

6 FORMAT! 1H0. 'BLANKET TENSION ='»F10.2,' LB/BLANKET ('.F8.1.' N/BL 
1ANKET) * ) 

7 FORMAT! 1H0» 'BUCKLING LOAD s»,F10.2»» LB ('»Fe.l*» N)») 

8 FORMATdHO. 'APPLIED BOOM LOADs. .F10.2. » OF BOOM BUCKLING LOAD') 

9 FORMAT (3F10.5.F11.5) 

11 FORMATdHO. 'REQUIRED POWER =»»F10.0.» WATTS') 

12 FORMAT ( 1H0. *DIA TO THICKNESS ='»F10.0.' IN/IN') 

13 FORMATdHO. 'BLANKET UNIT WT. ='.F10.3»» LB/FT**2 C.F8.2.' KG/M 
1 ** 2 ) 1 > 

17 FORMATdHO. 'LEADER LENGTH s»,Fl0.2.' FT C.F0.3,' M) • ) 

20 FORMATdHO. 'BLANKET AREA s*,F10.1»» FT**2 C.F8.2*' M**2 

1) •) 

21 FORMATdHO. 'ARRAY LENGTH ='»F10.2.« FT C.F8.3.' M) • ) 

22 FORMAT (1H0. 'ARRAY WIDTH s*,F10.2.» FT <*»FB.3,' M) • ) 

23 FORMATdHO. 'BOOM DIAMETER s»,F10.2»» IN C.F8.4.' M)«> 

24 FORMATdHO. 'BOOM THICKNESS s>,F10.3»» IN ('.F8.5»» M) ' ) 

25 FORMATdHO. 'BOOM EFFICIENCY s»,Fl0.2.' EFFECTIVE (I>») 

28 FORMATdHO. 'INITIAL BOOM (El ) = • rFiO.O . ' LB FT**2 ('.F8.0,' NM** 

12) •) 

29 FORMATdHO. 'PERCENT BUCKLING ='»F10.3> 

31 FORMATdHO. 'REQUIRED FREQ =»,F10.3»» HZ*) 

32 FORMATdHO. 'BOOM MATL DENSITY^' .F10.1. * LB/FT**3 C.F8.0.* KG/M 
1**3) ') 

33 FORMATdHO. 'BOOM MODULUS (E) ='.E10.4.' LB/FT**2 C.E0.3.* N/M* 
1*2)') 

35 FORMATdHO. 'BOOM UNIT WEIGHT ='.F10.2.» LB/FT <'»F8.2»* KG/M 

1) ») 

37 FORMATdHO. 'ACTUATOR WIDTH s»,F10.2.» FT C.F8.3.' M)») 

40 FORMAT(//1HO. 'CALCULATED SYM. FREQ.s* ,F7.3» • HZ.'.' CALCULATED A 
1SYM. FRE0.=*.F7.3.« HZ.') 

41 FORMATdHO. 'BLANKET WIDTH s*,F10.2.» FT C.F8.3.* M)«) 

42 FORMATdHO. 'BLANKET LENGTH S'»F10.2.» FT C.F8.3.* M) ' ) 

49 FORMAT ( 1H »'***♦* POWERs' ,F6 . 0 . ' WATTS AREA=» .F6. 1 . * FT**2 

1 WIDTH=' »F6. 3. * FT *****»/ih .»♦**** BOOM STIFFNESS=*. 
2F7.0» * LB*FT**2 TENSIONS' .F6. 2. • LB **♦**') 

50 FORMAT (1H '•**♦** POWERS' . F6. 0 . • WATTS AREAS' »F6.1 . * FT**2 

1 WIDTH=». F6. 3.' FT *****'/ih ,»**•** BOOM STIFFNESSs', 

2F7. 0. * LB*FT**2 *♦***♦) 

51 FORMAT ( 1H .'**♦** POWERS' ,F6.0» » WATTS AREAS », F6 . 1 » * FT**2 

1 WIDTHS' »F6. 3. ' FT *****'/iH .'*♦*** REQUIRED FREQUENC 

2Y=».F6.3.' HZ ***•*') 

249 FORMAT ( 1H »•**♦** POWERS' »F6. 0 . » WATTS AREAS '. F6 . 1 . ’ M**2 

1 WIDTHS* »F6. 3. ' M ***** '/1H , » ***** BOOM STIFFNESSs*. 

2F7 ■ 0 . ' N*M**2 TENSIONS* .F6. 2. * N ***♦*») 

250 FORMAT ( 1H .'****♦ POWERS' ,F6.0» * WATTS AREAS' , F6. 1 » * M**2 

1 WIDTHS* »F6. 3. ' M *****»/!H .'****♦ BOOM STIFFNESSs* . 

2F7 . 0 » ' N*M**2 ***♦*•) 

251 FORMAT ( 1H »»***♦* POWERS • »F6. 0 . * WATTS AREAS' ,F6. 1 » * M**2 

1 WIDTHS' »F6. 3. ' M *****'/lH »'♦**** REQUIRED FREQUENC 

2Ys* .F6*3. 1 HZ **♦**•> 

59 FORMAT ( 1H .'**♦** POWERS* ,F6. 0 > * WATTS BLANKET L='.F6.2»* FT 

1 BLANKET W='.F6.3»* FT *****'/lH .'****♦ BOOM STIFFNESSs', 

2F7.0 » * LB*FT**2 TENSIONS* »F6. 2. ' LB »****») 

54 FORMAT (1H .'**♦** POWERS' ,F6. 0 . • WATTS BLANKET L='»F6.2.» FT 

1 BLANKET WS*,F6.3»* FT *****«/ih ,»♦**** BOOM STIFFNESSs', 

2F7 • 0 » * LB*FT**2 *♦***•) 

61 FORMAT ( 1H .****** POWERS' ,F6. 0, » WATTS BLANKET Ls',F6.2»» FT 

1 BLANKET WS',F6.3»* FT *****'/ih ,»****• REQUIRED FREQUENC 

2Ys« ,F6*3. * HZ ****♦•) 

259 F0RMAT(1H »•**♦** POWERS* ,F6.0. ♦ WATTS BLANKET L='»F6.2.' M 

1 BLANKET W=',F6.3»* M ******/ih ,'**♦** BOOM STIFFNESSs'. 

2F7.0#* N*M**2 TENSIONS', F6. 2. * N •****») 

260 FORMAT ( 1H »•***♦* POWERS' ,F6.0, ' WATTS BLANKET L='»F6.2»' M 

1 BLANKET W=»,F6.3»» M *****'/ih ,»♦**** -BOOM STIFFNESSs'. 

2F7.0 . * N*M**2 ***♦*•) 

261 FORMAT ( 1H .'**♦** POWERS », F6 . 0 » * WATTS BLANKET L='»F6.2.' M 

1 BLANKET W=*,F6.3»* M *****»/ih ,»****♦ REQUIRED FREQUENC 

2Ys» , F6.3 » * HZ *♦***') 

52 FORMAT (////1H0 » ***** INPUT PARAMETERS FOR DEPLOYED ROLL-UP SOLAR AR 
1RAT****' ) 

53 FORMAT (1H1»10X,'****PRINTOUT OF PROGRAM OPTIMIZATION ITERATIONS*** 
1*»/22X» ' (ENGLISH UNITS-LB, FT * SEC .HZ) * ) 

56 FORMAT (1H » ******************************************************* 
1 **************************' ) 

57 FORMAT ( 1H1 . ' ****************************************************** 
X*************************** ) 

60 F0RMAT(1H1» » **********DEPLOYED ARRAY PARAMETERS**********') 

62 FORMAT!//* MODAL ANALYSIS FAILURE'/' TOTAL TENSION IS LARGER THAN 
1THE BOOM BUCKLING LOAD 0F'»F7.2,» POUNDS') 

63 FORMAT (//////1H »'*WARNING* TENSION IS BELOw MINIMUM RECOMMENDED T 

1ENSI0N 0F»/1H ,' 0.5 LBS/FT OF WIDTH REQUIREO TO ROLL UP 

2BLANKET ' ) 

NAMELIST/DATA1/NSETS .LEVEL. POWER . AREA . EIBOM1 .ROARAY . WIOTH.EIEFF. TL 
1EADR.FREQR.RBOOM. DPT .EBOOM. TENS. NDATA2.NUNITS.BKTW.BKTL.NTYPE.PCTB 
2KL 
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c 


ALL VARIABLES STARTING WITH THE LETTERS GE ARE PROTOTYPE ARRAY VALUES 
OATA GEEIBE,GEROBE»DEGOP*GEDlAB.GEWlDC,/3040.» .11*20. *1.34» .5/ 

1 GE w» GET/8.25 »2./ 

MAXITR=15 
NW-2 
NL=5 
LOP-0 

100 LOP=LOP+l 
READ (5»DATA1 ) 

RPT=DPT 

BKLPCT=PCTBKL 

PCTBKL=ABS(PCTBKL> 

IF INTYPE. EO.O)AREA=BKTL*BKTW*2. 

IF(LEVEL.EQ.l)GO TO 111 
IF(LEVEL.EQ.2>G0 TO 111 

IF(LEVEL.EQ.3)G0 to 111 
IF(NTYP£.NE.O)60 TO 102 
IFtBKLPCT.LT. 0.0)60 TO 101 

TFNS = 3.l4l59**?*FIBOMl*t>CT?KL/( 20 rt .*(B*TL+Tl FA OR ) **? ) 

GO TO 111 

101 E I 50“ 1 *200 , *TEN5* ( BKTL+TLEA09 )**2/ ( 3 . i Ai?o**p*p<-TB<L ) 

GO TO 11! 

1 02 TF(NUNITS.FO.O»GO TO 104 
Pf ROMlsETBOvil/.4l325331 
EBOOM*EBOOM/47. 880259 
W!DTH*WI0TH/.3C48 
TENS a TF.NS/4 ,4482216 
AREA* AREA/. 09 290304 
TLEADRsTLEADR/. 3048- 

104 IF(BKLPCT.LT.O.O) GO TO 103 

0! ABOMs (2»*(FIB0Ml*RPT/( FROOM*F |FFF*4.*3.l4l59))**.25)*12. 

WI0CEs(DIAE0M/GEDIAB1*GFWI0C 

L*AREA/ (WlDTH-WIDCEl+TLFADR 

TFNS*3« 141 59**2 *R!90Ml*PCTBKL/{ 20 O.*l**2) 

GO TO 105 

1 02 TALL RSAZR01F1B0MI .TFNS.2) 

105 IFfNUNFTS.EO.OlGO TO 111 
AREA= AREA*. 09290304 
TLEADR=TLEADR*.3048 
TENS=TEN5*4. 44822 16 

Wl 0TH=W1DTH*. 3048 
E“OOM=EBCOM*47. 880259 
F I 90M1=E I B0M1 *.41325331 
111 TF(NTYPE.EO.O)GO TO 108 
WR 1 TE ( 6 » R 7 ) 

WRITF16.56) 

IF (NUN I TS.EO. 1 )60 TO 106 

IFtLEVEL.EO.l > WR I TE ( 6. 51 I POWER. AREA. WI nTH.FRFOR 
IF(LFVFL.E0.2 1WRI TE( 6. 50) POWER, AREA, w] f>TH,FlPOMl 
I F f IABStLEVEL).E0.3)WRI TE( 6.49) POWER .AREA. W I OTH. El 30M 1 , TENS 
GO TO 107 

106 IFCLEVEL.EO.l >WRI TE(6»251 ) POWER .AREA .WIOTH.FREOP 

l F ( LEVEL. EQ. 2 )WRITE(6»250) POWER « AREA , w l OTH, F l R0M1 
I F ( I ABS(LF.VEL). F-0. 3 ) WRITE ( 6.249 ) POWER, AREA, WIOTH.EIBOMI .TENS 

107 WRITE(6,56) 

WRITE (6,56) 

GO TO 114 

108 wRItE(6,57) 

wR I TF ( 6 , 56 ) 

IF (NUN! TS.EO. 1 )G0 TO 109 

IF(LEVEL.EQ.l )WRITF(6,61 > POWER, 9KT1., BKTW.FPFQR 
I F ( LEVEL. E0.2 ) WRITE (6, 54) POWER, 0KTI.,B5TW,E!BOM1 
I F ( I ABS( LEVEL ).E0. 3) WRITE (6,59) POWER.BKTL.BKTW, El BOM 1 .TENS 
GO TO 110 

109 IF (LEVEL. EQ.l) WRITE (6 » 261) POWER *BKTL»BKTW.FREQR 
IF (LEVEL. EQ. 2) WRITE (6. 260) POWER .BKTL'BKTW.EIBOMl 
IF(IABS(LEVEL) .EO. 3) WRITE (6, 259) POWER »BKTL*BKTW»EIBOMlr TENS 

110 WRITE ( 6’ 56) 

WRITE ( 6 » 56) 

114 IF(NUNITS.EQ.0)GO TO 115 

c 

C SINCE ALL CALCULATIONS ARE PERFORMED IN ENGLISH UNITS* 

c this section converts units if the user inputs si units, 
c 

ARFA=AREA/. 09290304 
TLEADR=TLE ADR/. 3048 
TENS=TENS/4. 4482216 
ROAR AY=R0ARAY/4. 88242764 
• WIDTH=WI0TH/. 3048 
R<TL=BKTL/.3048 
RKTW=BKTW/. 3048 
RPOOM=RBOOM/ 16. 0184637 
FROOM = EBOOM/47. 880259 
EIB0M1=EIB0M1 /.41 325331 

115 WRITE16.52) 

WRITE (6,11 )POWER 
SIUNIT=AREA*.092903C4 
WRI TE< 6,20 ) AREA, SI UN IT 
SIUNIT=WIDTH*.3048 

I F INTYPE. EO.l )WR I TE( 6*22) WIDTH, SlllNI T 
«IUNlT=BKTW*.3048 

IF(NTYPE.E0.n)wRITE(6.41 )BKTW.SIUNIT 
SI UNIT=BKTL*. 3048 

JF(NTYPE.EO,0) WR I TE ( 6 » 42 ) BK TL.SI UN IT 

SIUNIT=TLEADR*. 3048 

WRITE(6,17)TLEADR.SIUNIT 

S I UN I T = R0ARAY*4 . 88242764 

WRITE (6,13)R0ARAY.SIUNIT 

I F < JABS (LEVEL ).LE.l ) WRITE (6.31) FREQR 

S I UN I T = E I BOM 1*. 41325331 

I F ( LEVFL. EO. 1 ) WRITE (6. 28 )EI80M1 .Si UN I T 

IF( IA3SI LEVEL ).GF. 2) WRITE (6.4) El B0M1 • STUN IT 

SI UN I T a E ROOM* 4 7 ,880259 

WRITE ( 6 , 33 ) E BOOM , S 1 UN I T 

« I UN I T = RROOM* 16.0184637 

WRITE (6,32 )RBOOM*SI UNI T 

WRITE <6,25)£IEFF 
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WRITE ( 6* 12 )RPT 

IF(LEVEL.EQ.-3)WRITE(6*29)PCTBKL 
I F { I A6S( LEVEL ).NE.3 )G0 TO 205 
EIROOM=EIBOM1 

01 ABOM=(2.*(EIBOOM*RPT/(EBOOM*EI£FF*4.*3. 141 59 > >**.25) *12. 
THIB0M=DI ABOM/RPT 

ROBOOM=4.*(3. 14159*EI800M/(E800M*EIEFF*RPT ) 1 **• 5D*RB00M*( • -DFGO 
IP) /360. 

WI DOE* ( DI A BOM /G ED I A B ) *GEW IOC 

I F (NTYPE.EQ.l »WPANEL=W1 DTH-WIOC E 

IF(NTYPE.E0.0»WPANEL=B<TW*2. 

!F(NTYPE.EO.O)WIDTH=wPANEL+WIDCE 
I F( NTYPE.EQ.l )L=AREA/WPANEL+TLEADR 
IF,(NTYPE.EQ.O )L=BKTL+TLEADR 
SI UNIT =TEN5*4.44822l6 
WRITE <6 *6) TENS. SI UNIT 
T 1=3 • 1 4159**2*E I80M1 /L**2 
IFlTFNS.LT. 0.5*T1 ) GOTO 150 
WRITEI6.62) T 1 
GOTO 999 

150 FDIF=RSAFX(TENS»EIB00M*FREQA,FREQS) 

I F ( FREQS.EO.O • ) GOTO 999 
GO TO 211 

205 WRITE(6.*3) 

IF (LEVEL.NE.2 1G0 TO 215 

E I BOOM = E 1 B0M1 

CALL RSAZRF(TENS.EIBOOM) 

GO TO 211 

215 CALL RSAZER(EIB00M.EIB0M1,TENS> 

01 ABOM=(2.*(EIBOOM*RPT/(EBOOM*EIEFF*4.*3.14l59) )**.25)*12. 
THIBOM=DIABOM/RPT 

R0p00M=4.*{ 3.14159*EIB00M/(E800M*EIEFF*RPT ) ) ** . 5O*RR0OM* ( 360 . -DEGO 
IP) /360. 

WI DCE= ( DI ABOM/GEOI AB) *GEW! DC 
IF (NTYPE.EQ.l )WPANcL=WlDTH-WlDCE 
IF< NTYPE. EQ.0)WPANEL=BKTW*2. 

I F(NTYPE.EO.O)WIDTH=WPANEL+W I DC E 
! F (NTYPE.EQ.l >L*AREA/WPANEL+TLE ADR 
IF(NTYPE.EQ.O)L*BK,TL+TLEADR 

211 WR I TE ( 6 » 60 ) 

WRITE (6*11 ) POWER 

SI UNIT= ARE A*. 09290304 

WRITE (6. 20) AREA. SI UN IT 

I F (NTYPE .EO. 1 )BKTL=AREA/WPANEL 

T F ( NTYPE • EQ. 1 )BKTW=WPANFL/2. 

Si UN I T=BYTW*. 3048 

WRITE(6.4l)BKTW»StUNIT 

S1UNIT=BKTL«.304B 

WR I TE < 6 *42 ) BKTL *S IUNI T 

SI UN I T=WIDTH*.3048 

WRITE (6.22)WIDTH.SIUNIT 

SIUNIT=L* .3048 

WRITE (6.21JL.SIUNIT 

S I UN I T*W I OCE* . 3 048 

WRITE(6.37)WI0CE.SIUNIT 

SIUNIT=TLEADR*.3048 

WR I TE( 6 . 1 7 ) TLEADR.SI UNI T 

S I UN I T = ROAR AY *4 . 8824 2764 

WRITE (6*13) ROARA Y.S IUN I T 

SIUNIT=TENS*4.4482216 

WRITE(6,6)TENS.SIUNIT 

BUCK=3. 14 1 59**2*E 1 BOOM/L**2 

BOOMl D=2.*TENS/BUCK 

WRITE(6»8) BOOML 0 

S I UN I T = BUCK *4 .4482216 

WRITE(6,7) BUCK.SIUNIT 

STUN I T = EI BOOM*. 41 32 5331 

WRITE < 6.4 JEIBOOM.SIUNIT 

SI UN I T = DI A BOM* • 0254 

WRITE ( 6.23)01 AB0M.SIUN1T 

SI UNIT=THI BOH*. 0254 

WRITE ( 6 *24 )T HI BOM* SI UN IT 

SIUNIT=R0B00M*1 .4881639 

WRITE(6»35) ROBOOM.S I UNIT 

S I UN I T =EB00M*47 .880259 

WRITE (6.33 ) EBOOM » S IUN I T 

S I UN I T = RROOM* 16 .018463 7 

WRITE (6*32 )R BOOM, 5 1 UNIT 

WRITE (6»25)EIEFF 

IF ( I ABS( LEVEL ) • LF.2 ) WR I TE ( 6*5)FRE0 

I F ( I ABS(LEVFL) .GE.3 )WRITF(6,40) FREQS > r RFOA 

I F ( 2 . *TFNS/WPANEL .L T .0.5) WR I TE ( 6.63) 

F I B0M1 = E I BOOM 

IF ( NDATA2.NE.0 )CALL RSABAS( TENS *E I BOOM) 

PCTBKL=BKLPCT 
I F ( NUN I TS • EQ. 0 ) GO TO 998 
AREA = AREA*. 09290304 
TLEA0R=TLEADR*.3048 
TENS=TENS*4. 44822 16 
R0aRAY=R0ARAY*4. 88242764 
WI DTH=WIDTH*.3048 
oKtl=3KTL*.3048 
•RKTW=BKTW*.3048 
R BOOM =R BOOM* 16.0184637 
EB00M=FB00M*47. 880259 
F I R0M1 = FI BOM 1*. 41 325331 

998 I F ( LOP.LT .NSETS )G0 TO 100 

999 STOP 
END 
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SUBROUTINE RSAZER Ut *X1 *T I 

REAL X.Xl.EpS.AER. RE« » F» OF* OLQF. OLDDF.T OL* O* DX* OLOX »XX * A* B» C* FA# FB 
1 »FC .HE AN »L 

COMMON/ JSP1/ EIEf F. FR EQ. 01 ABQN.T HIBOH, ROBOOM.W I0CE.0E60P. WPANEU MAX 
HTR.HDAT2.R0ARAf • NW *NL »L *W *F RE QR • RB OOK.EB OOH. RPT. POWER. AREA. TLEADR 
2.BKTW.BKTL.NTYPE 

11 FORMAT (/1HQ.*B00M STIFFNESS OPTIMIZATION ITERATIONS* /* I TE 

INS ION BOOM I £11 FREQ OIF FREQ FREQ REQM 

12 FORMAT (1H .12 .£E 14 .5 • IK* I 

32 FORM AT ( 45X. * BOOM STIFFNESS OPTIMIZATION DIOKT C0WER8EM 

IT£R=0 

EPS=FR£QR*.02 
AER-X1 • • 005 
RERZ.OIO 
X= XI 

CALL RSAZRFfT.X) 

F=FREQ-FREQR 
FREQlzFREQ 
WRITE C 6 .11 1 

WRITEIS.121ITER.T. X. F. FREQI. FREQR 

IFIABStF )*LE.EPSI RETURN 

IT£R=1 

TOL=AER/RER 

D=( ABS(X.RER) ♦AER) *5 .0 
CALL RSAZRF (T.X*D I 
DF=l IFREQ-FREQRI-Fl/D 
60 TO 20 
IS WRITE(6.11I 

WRITE (6.12 IITER.T »X .F . FREQI .FREQR 
18 ITER = ITEft«-l 

D=(ABS (X«RER)«AER 1*5.0 
CALL RS AZRFCT.X*D) 

DFr | CFREQ-FREQR l-F) /D 

IF( OF*OL DDF. GT.O •) 00 TO 20 

OX=X-OLDX 

0 = 2 . 

GO TO 25 
20 DX=-F/DF 

XX=DX/ (SIGN (TOL »X »♦ X I 
0=AM AX1 ( 0 « 5 • XX. - 2. *X X I 
2 S OLD X= X 
OLOF=F 
0LDFRQ=FREQ1 
OLOOF=DF 
X=OLDX*f)X/D 
CALL RSAZRF(T.X) 

FZFREQ-FREQR 

FREQ1=FREQ 

IFIABSCF 1.LE.EPS) GO TO 85 
IFIF • OLDF.LE.O.) ,G0 TO 34 
IF I ITE R.LT.MAXITR > GO TO 15 
PRINT 32 
ITER=-MAXITR 
GO TO 85 

34 IF ( ABS (DX/D 1.LE.ABS (X«RER>+AER) 60 TO 05 

A=Q10X 
B = X 

FA=OLDF 

FB=F 

AFREQ=0LOFRQ 
BFREQ=FREQ1 
60 TO 60 

40 WRITE (6.11 ) 

WRITE(6.12)ITER«T. B. F8. BFREQ. FREQR 
IF(ABStX-fi).LT.TOL) X =B* SIGN ( TOL *C-B I 
A=B 
FA = FB 

AFREO=BFR EG 

IF( (X-MEAN )*(8-X) .GT.O.I GO TO 50 

b=mean 

GO TO 55 
50 B=X 

55 CALL RSAZRF (T.B 1 
FB=FREQ-FREQR 
6FREQ=FREQ 
ITER = n ER*1 

IF(FC*F8.LT.O.) 60 TO 65 
60 C=A 

FC=F A 

CFREQZAFRE6 

65 IF( ABSI FB J «L T • ABSf FC 1 1 60 TO 70 

A=B 
FA=FB 

AFRE0=8FREQ 

B=C 

FB=FC 

BFREQZCFREQ 

C=A 

FC=F A 

CFREQ=AFRE 0 

70 IFf A8S( FBI »L T.EP SI GO TO 80 
MEAN =(B*C)/2. 

IF( FB.N E.FAI X = ( A*FB- B* FAJ / (FB- FAI 
IF ( FB » EO.FA I X = It AN 
TOL = ABSCB»RER»*AER 
IFC ABSCNEAN-8) .LT.TOLI 60 TO 80 
IF( ITER «L T.M AXITR) GO TO 40 
PRINT 32 
IT ER=-M AXITR 
0 0 X=B 
F=F8 

FREQZBFREQ 
85 WRlTEt6.ll) 

WRITE (6.12) ITER.T.X.F .FREQ. FREQR 

RETURN 

END 
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SUBROUTINE RSAZRFt X* EIBOOM) 

REAL X #X1 »EPS»F *A #B *C » FA #FB »FC »HEA N »L 

COHMON/JSP1/EIEFF, FREQ. DI ABOM# T HI BOM# R0800M#y IDCE# DEGOP# UPANEL. H AX 
1 ITR *ND AT2 » ROAR A Y » NU #NL *L*W #F RE QR t RB OOK.EB OOH» RPT • POWER* AREA ■ TLEAOR 
2w BKTU.BKTL.NTYPE 

COMMON AJSP2/GEDIAB# GEN IDC#GEW» GET. GEEIBE.GE ROBE 

11 FORM ATI//SX. 'TENSIOH 0? T IMIZ ATION IT ERAT IONS * /5X * I TENSION 

1 BOOM (ED FREQ DIF SYM FREQ ASTH FREOM 

12 FORMATtlH »SX*X2.SE14.S> 

24 FORMAT I50X » * END TENSION OPTIMIZATION* I 

32 FORM ATI sox# * TENSION optimization oidnt C0NVER6E*) 

39 FORMAT I50X»* NO ROOT IN BOUND* 1 
IFI EIBOOM .LE.O.) 63 TO 100 

DIABOH = I2.*IEIBOOK*RPT/IEBOOH*EIEFF*4.*3.14153 I)**. 25 1*12. 
THI0OM=DI ABOM/RPT 

R0B00M=4.* <3.1415 S*EIB OOH/(EBOOM«EIEFF *RPT 1 1 ••.50 • RBOOH*f 360. -0E60 
1PI/360. 

W IDCE- IDIABOM/GED IAB)*GEWIDC 
IFINTYPE.EQ.DUP an el=u-widce 
ifintype.eq.0)upanel=bktu*2. 

IFINTYP E. EG.O)U=U? A* EL*W IDCE 

IF< NTYPE .£ 0.1 ILSA PE A/W PANE L+TLE ADR 

IFINTYP E. EQ.0>L=B<TL*TLEA0R 

Xl=l .0 2*3.14159*»2*£IB00H/L**2>/2. 

X2=f .98*3.141£9* *2 •EIB00M/L** 2) / 2. 

EPS=.005*FREQR 
ITER=1 
A -XI 
8=X2 

FA = RSAFX CAtEIBOOM »F R£QA»FR£GS) 

AFREQA=fR EGA 
AFREQS=FREGS 

F8=RS AFX I B. EIBOOM# FR EQA.FREQS) 

BFREQA=FRE GA 
BFREQS=FREQS 

IFIFB*FA.GT.0.1 GO TO 36 
MRITEI6.11I 
60 TO 60 

36 PRINT 39 

GO TO 101 

40 URITE<6.12)ITER.8» EIBOOM# FB. BFREQS. BFREQ A 
45 IFIABSIX-B J.LT.TOL) X =B ♦ SIGN I TOL .C-B ) 

A=8 

FA=F8 

AFREQ A= 8FREQA 
AFREQS=BFREQS 

IFUX-N£AN1*<B-X).GT.0.) 60 TO 50 
8-NFA N 
GO TO 55 

50 BIX 

55 FB=RSAFXlB#EIBOOM. FREQA.FREOSt 
BFREQA=FRE QA 
BFREQS=FREQS 
ITER=ITER»1 

IF< FC*FB-LT.O.) G3 TO 65 
60 C-A 

FC-FA 

CFREQA =AFREQA 
CFREQS = AFREQS 

65 IFIABSIFB).LT.ABS«FC)> GO TO 70 
A- B 
FA = FB 

AFREQAr BFREQ A 

AFREQS=BFREOS 

B=C 

FB=FC 

BFREQ A= CFREQA 

BFREGS=CFREQS 

C=A 

FC-FA 

CFREQ A= AFREQ A 
CFREQS=AFREQS 

70 IFI ABSC FB) .LT.EPSI GO TO 80 
MEAN - f B *C 1/2. 

IFIF8.NE.FA) X = ( A« FB- 8* FA) / IF B- F A) 

IFIFB.EQ.FA I X = ICAN 
TOL = ABS IB*. 001) 

IFIABStKEAN-BI .LT.TOLI GO TO 60 
IFI ITER .LT.MAXITRI GO TO 40 
101 PRINT 32 

IT ER--M AX ITR 
8 0 X=B 

F-FB 

FREQS=BFREQS 
FREQ A=BFR EGA 
FREQ=(FR£GS*FREQA 1/2. 

UR IT El 6* 12 ) ITER# X# EIBOOM# F. FREQS. FREQ A 
GO TO 26 

100 FREQ = El BOOM 
2 6 URI TE I 6 #24 ) 

RETURN 

EHD 
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FUNCTION RSAFX1T.E IBOON. FREQ A. FREQS » 

REAL H(4Q00).K 14000 >.VALC100).L 

COMMON/ JSP1/ EIEFFt FR EQ. DIA80M.THIB0N, ROBOOH.V IDCE.DE60P. UPANEU WAX 
1 1 TR > NO AT 2 • ROAR A Y » MU » NL #L * W .F RE QR • RB OOM »EB OOK.RPT. POWER. AREA* TLEADR 
2. BKT V. 8KTL *NTYP E 

COMMON /JSP 2/GEDIA B* SEWIOC.GEW.GET »6EE IBE * GE ROBE 
COMMON/ JSP6/VL A* AC • V L AMT • VL AMH* RHOBEt AKBE.S IGMBE » EBE* VLMOIA 
3 FORMAT** MODAL ANALYSIS FAILURE*) 

21 FORMAT* 2X .1H»*IS) 

W2=W/2. 

yp ANE2 = WP AN El/2. 

EIBEAN=6EEIBE* (O. «• T/OET^O .5 *V LA MAC )«W/6EW 
R0BEAN=6ER0 9E*SQRT IE IBE AN/ GEEI8E) 

T2=-2.*T 

RO AR AY-RO AR AY • I MP ANE2/y 2 ) 

CALL KRAYIK»N»EIBEANtEIBOOM«T»T*T2*TZ»NU«NL»L»U2tO ) 

CALL MR AY IH » N *R0 BE AM »R 0 BOOM* RO AR AY « NW • NL* L.y 2 *0 ) 

CALL GENIGNIK.K.VAL *VA L»-N.O.OKI 
ROAR AY=R0ARAY*(U2/ y? ANEZ) 

IF* OK.EQ.O. »G0 TO 1 
DO 10 1=1. N 
A = VA LII) 

IF * A.GT.O.) GO H 15 

10 CONTINUE 

15 IF* I.6T .1 »yRITE<6. 211 1 
FREflS= SORT *A)/6. 283185 3 
T=T* ( U2 / y P AN £2 ) 

T2=-2. «T 

CALL K.t ay (X.N.EIBE am . El BOOM. T.TtT2.T2.Ny.NL. L.y 2.1) 

CALL MRAYfM.N.ROBEAN.R OBOOM.ROARA Y.NW.NL.L.W2.1 ) 

CALL GENlGN (M.K.VAL.VAL.-'N.O.OK) 

T=T • (U PANE 2 /y 2 ) 

T 2 =- 2 . • T 

IF* OK.EQ.O. IGO TO 1 
00 11 1=1 .N 
A=VAL ( I) 

IF * A.GT.O.) GO T) 16 

11 CONTINUE 

16 IFlI.GT .4)URITE16» 21 I I 
FREQA= SORT (A )/6 .28318 5 3 
RS AFX=FREOS- FREQ A 

GO TO 2 

1 RS AFX=T *100. 

F RE 0S= 0. 

WRIT EI6.3 t 

2 RETURN 
ENO 


SUBROUTINE RSA8AS ( T * E I BOOM ) 

REAL L 

C0MM0N/JSP1/EIEFF, FREQ, DIABOM,THIBOM.ROBOOM,WIDCE»DeGOP»WPANEL. MAX 
1ITR, NDAT2,ROARAY,NW,NL,L#W*FREQR,RBOOM,£BOOM,RPT, POWER, AREA, TLEADR 
2»BKTW.BKTL»NTYPE 

COMMON/ JSP2/GEDI AB . GEW IDC , GEW , GET , GEE IBE , GEROBE 
COMMON/ JSP5/C0NST . COEFF , EXP 

COMMON/ JSP6/VLAMAC , VLAMT . VL AMH » RHOBE » AKBE , SI GMBE . EBE , VLMDI A 
1 FORMAT*/* TOTAL BLANKET WEIGHT =*,F9.2»* LBS *»,F7.2»' KG) 

1 *) 


2 

FORMAT* *0 

DRUM SHELL WT, 

=*,F7.2,* 

LBS 1* 

#F6.2» 

* KG) * ) 


3 

FORMAT* • 

END CAP+GUIOE WT. 

=*,F7.2,* 

LBS (• 

, F6.2 » 

' KG)*) 


4 

FORMAT** 

SUPPORT SHAFT WT. 

=*»F7.2»* 

LBS *» 

•F6.2* 

• K6) • ) 


5 

FORMAT* * 

BEARING ASSY WT. 

=* ,F7 .2, * 

LBS *• 

»F6.2* 

• KG ) * ) 


6 

FORMAT*/* 

1»> 

FORMAT*/* 

LEADING EDGE BEAM WEIGHT = 

»»F9.2* 

• LBS 

(• ,F7.2* ' 

KG) 

7 

OUTBOARD END SUPPORT 

ASSY WT.= 

• ,F9.2# 

» LBS 

C * »F7 .2 # ' 

KG) 

8 

FORMAT* *0 

CENTER SUPPORT WT, 

=» ,F7.2»* 

LBS *♦ 

»F6.2» 

• KG) * ) 


9 

FORMAT* » 

SLIPRING+HARNESS WT=»»F7.2,* 

LBS 

»F6.2, 

• KG)*) 


10 

FORMAT* /• 
1*) 

FORMAT* » 

DEPLOYED BOOM WEIGHT 

= 

*,F9.2»* LBS 

f «,F7.2»* 

KG) 

11 

NEGATOR ASSY WT. 

=*,F7.2,* 

LBS ** ,F6.2» 

• KG)*) 


12 

FORMAT* » 

BOOM ACTUATOR WT. 

=*.F7.2»* 

LBS 1* 

»F6.2» 

« KG)*) 


13 

FORMAT*/* 

1»> 

FORMAT* » 

TOTAL SOLAR ARRAY WEIGHT = 

* ,F9.2, 

* LBS 

(• »F7.2» * 

KG) 

14 

TOTAL ORUM ASSY WEIGHT = 

• , F9.2* 

» LBS 

* * ,F7.2» * 

KG) 

15 

FORMAT* • 

total CENTER support 

ASSY WT = 

».F9.2» 

» LBS 

* * ,F7.2» * 

KG) 


1 *) 

41 FORMAT (2110) 

42 FORMAT < 9F8.4 ) 

43 FORMAT*/* POWER TO WEIGHT EFFICIENCY =',F9.2»* WATTS/LB * • »F7.2» 
1* WATTS/KG)*) 

47 FORMAT * ///1H0 » » *******C ALCULATED BASE STRUCTURE PARAMETERS******* 

1*/* (LAUNCH ACCELERATION=* ,F5.2r **BASELINE LEVEL)*) 

50 FORMAT (1H1»*****BASE STRUCTURE INPUT PARAMETERS*****) 


100 

FORMAT *10F7 

3) 





101 

FORMAT *1H0» 

RHOS 

ES 

SIGMS RHOC 

EC 

SIGkC*) 

102 

FORMAT *1H0» 

RHOSS 

ESS 

SIGMSS RHOBE 

EBE 

SI GMBE*) 

103 

FORMAT *1H0# 

RHOES 

EES 

SIGMES RHOCS 

ECS 

SIGMCS * ) 

104 

FORMAT *1H0# 

VLAMAC*) 






NAMELIST/DATA2/RH0S.ES»SIGMS»EC»SIGMC.RH0C*RH0SS»ESS»SI6MSS»EBE*SI 
IGMBE » RHOBE ♦ SI GMES , EES , RHOES , S I GMCS , RHOCS » ECS * VLAMAC 
DATA RH0S»ES»SIGMS,AKS»EC#AKC»SIGMC»RH0C/1..1.*1.,1.»1.#1. »1.»1./ 
1RH0SS»ESS»AKSS,SIGMSS#RH0B#EBE*AKBE/1.*1.»1.»1.»1.#1.»1./ 

2SIGM8E » RHOBE# AKES.SIGMES, EES# RHOES, AKCS»/1.,1. #1. »1.»1.#1. »1./ 
3SIGMCS , RHOCS # ECS • RHOA » AK A # E A * SI6MA , RHONEG/1 .*l.*l.,l.,l.*l.»l..l./ 
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4VLAMAC/1./ 

DATA GEWS»GEWC»GEWSS»GEWB»GEWBE»GEWES/5.7B»4.53»1 .21*2.50* .85*4.1/ 
lGEWCS»GEWNSfGEWA*GEWNEG.6EBLKT/1.33*2.34»5.ie*i.54»46.50/ 

2GEB00M » GETOTL # GER I » GERO » GEL/6 . 55 * 82 . 5 » 4 . » 5 . * 33 . 2/ 

acc=vlamac 

NBAS=0 

6EPANE=6EW-GEWIDC 
97 READ(5.DATA2) 

ACC—Vl AMAT 

99 VLMDIA=DIAB0M/GEDIAB 
VLAMT=T/GET 
VLAMH=WPANEL/GEPANE 
VLAMV=L/GEL 
VLAMHP=W/GE* 

WRITE (6»50 ) 

WRITE (6*101) 

WRITE (6t 100) RHOS’EStSlGMSi RHOCtECrSIGMC 
WRITE (6* 102) 

WR I TE < 6 » 1 0 0 ) RHOSS * ESS » SI GMSS * RHOBE* EBE » SI GMBE 
WRITE(6» 103) 

WRITE (6» 100) RHOES » EES » SIGMES » RHOCS » ECS » S I GMCS 
WRITE (6* 104) 

WRITE(6»100)ACC 
WBOOM=ROBOQM*L 
WBLNKT=ROARAY*L*WPANEL 
wns=gewns*vlamhp*vlamv 
WNEG=GEWNEG*RHONEG*SGRT ( VLAMT ) 

RATIO=SQRT(l./(VLAMH*VLAMV)) 

300 VLAMAC=ACC*RATIO 
WST=GEBLKT*GEWS 

const=wblnkt/wst 

RATIOW-GEWS/WST 

COEFF=RATIOW*<RHOS*VLAMH**(5./3.>*VLAMAC**(2./3.)*SQRT<1./(SIGMS*A 
IKS) )/ES**(l,/6. > ) 

EXP=2./3. 

CALL RSAZRO(WBSR»T#0> 

VSST= ( VLAMAC*VLAMH*WBSR/ES ) ** <1 . /3. > 

VSS0=(VLAMAC*VLAMH*WBSR)**(l./3. )*SQRT(ES**(l./3.)/(SIGMS*AKS) ) 
WS=WBSR*WST-WBLNKT 

VCST=SQRT( (VLAMAC*WBSR/<AKC*SIGMC> ) ) 

WC=GEWC*RH0C*VSS0*VCST**2 

wst=geblkt+gews+gewc+gewss 

CONST=( (GE8LKT+GEWS)/WST>*WBSR-fRH0C*VSS0*VCST**2»GEWC/WST 
RATIOW=GEWSS/WST 

COEFF=RAT I OW* ( RHOSS*VLAMH*» < 5./3. ) * VLAMAC** ( 2 . /3 . ) *SQRT ( 1 . / (SIGMSS 
1+AKSS) >/ESS**(l./6.) > 

EXP=2./3. 

CALL RSAZRO ( WBSCSR » T • 0 ) 

VSSST= t VLAMAC*VLAMH*WBSCSR/ESS ) ** ( 1 . /3 . > 

VSSSO= ( VLAMAC* VLAMH*WBSCSR > ** ( 1 . /3 . > *SGRT ( ESS** < 1 • /3 . ) / ( SIGMSS* AKS 
IS)) 

wss=wbscsr**st-wblnkt-ws-wc 
WB=GEWB*RH 0B*VSSS0**2 
CALL RSAZRO(WBER»T»l> 

wbe=wber*gewbe 

VBEST= ( ( WBER*VLAMAC* . 5+VLAMT* . 5 ) *VLAMHP/EBE ) **( 1./3. ) 
VBES0=SGRT{(WBER*VLAMAC*.5+VLAMT*.5)*VLAMHP/(VBEST*AKB_E*SIGMBE)) 
VLMVES=SQRT( VLAMV+ < VSSO**2-VLAMV ) *GERI**2/GERO**2) 

VESST= ( VLAMAC* VLMVES*WBSCSR ) ** ( 1 ./3 . ) / < AKES*SIGMES*EES ) ** <1 . /6 . > 
VESS0=(VLAMAC*VLMVES*WBSCSR/(VESST**2*AKES*SIGMES) ) 
WES=GEWES*RHOES*VESST*VESSO*VLMVES 
WST=GEWA+GEBOOM 

const=wboom/wst 

RATIOW=GEWA/WST 

COEFF=RATIOW*RHOA*VLMDI A*VLAMV** ( 1 . /2 . ) * < VLAMAC/ ( AKA*SIGMA ) ) ** ( 1 . / 

12.1 

EXP=l./2. 

CALL RSAZRO(W8AR*T»0) 

WA=WBAR*WST-WBOOM 

VAST= C VLAMAC*WBAR/ < AKA*SIGMA ) ) ** ( 1 . /2 . ) 

WST=GEBLKT+GEWS+GEWC+GEWSS4-GEWA+GEBOOM+GEWB*-GEWNEG+GEWNS 
WBSCSA= ( WBLNKT+WS+WC+WSS+W A+W800M+WB+WNE6+WNS ) /WST 
VCSST=SQRT(VLAMAC*VLMVES*WBSCSA/(VLMDIA*AKCS*SIGMCS) ) 
WCS=GEWCS*RHOCS*VLMDIA*VLMVES*VCSST 
WT=W8LNKT+WS+WC+WSS+WB+WBE+WES+WCS+WNS*WBOOM+WA4-WNEG 
WDRM=WS+WC 

WBUS=WCS-*-WA*WB+WNEG+WNS+WSS 

WATPS=POWER/WT 

IF(ABS(RATlO-SQRT(GETOTL/WT)).LT.O.i*SGRT(GETOTL/WT))GO TO 301 
RATIO=SQRT ( GETOTL/WT ) 

GO TO 300 

301 WRITE (6* 47) VLAMAC 
SIUNIT=.454*WBLNKT 
WRITE(6*1)WBLNKT»SIUNIT 
SIUNIT=.454*WB00M 
WRITE (6» 10) NBOOM'SIUNIT 
SIUNIT=.454*WBE 

WRITE (b#6)WBE»SIUNIT 
SIUNIT=.454*WES 
WRITE (6»7 ) WES iSIUNIT 
SIUNIT=.454*WS 
WRITE <6*2)WS*SIUNIT 
SIUNIT=.454*WC 
WRITE<6»3) WC*SIUNIT 
SIUNlTs. 454* WORM 
WRITE(6*14) WDRM»SIUNlT 
SIUNIT=.454*WCS 
WRITE (6* B) WCS» SlUNIT 
SIUNIT=.454*WA 
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WRITE(6»12)WA»SIUNIT 

SIUNIT=.454*WB 

WRITE(6»5)WB»SIUNIT 

SIUNIT=.454*WNEG 

WRITE (6* 11) WNEG.SIUNIT 

SIUNIT=.454*WNS 

WRITE(6»9)WNS»SIUNIT 

S1UNIT=.454*WSS 

WRITE (6»4)WSS»SIUN IT 

SIUNIT=.454*WBUS 

WRITE (6t 15) WBUS» SIUNIT 

SIUNIT=.454*WT 

WRITE(6»13)WT*SIUNIT 

SIUNlT=WATPS/.454 

WRITE<6» 43) WATPSr SIUNIT 

199 NBAS=NBAS+1 
VUAMAC=ACC 

IF(NDAT2.GT.NBAS)G0 TO 97 

200 RETURN 
END 


SUBROUTINE RSAZRO ( X » T * NCHECK ) 

REAL MEANiL 

COMMON/JSPl/EIEFF » FREG » Dl ABOM» THIBOM • ROBOOM. WIDCE»DEGOP»WPANEL»MAX 
1 1 TR # NOAT 2 * ROAR AY » NW . NL . L # W » FREQR # RBOOM » EBOOM * RPT , POWER » AREA > TLE AOR 
2»BKTW#BKTL#NTTPE 
COMMON/ JSP5/CONST.COEFF»EXP 

COMMON/ JSP6/VLAMAC * VLAMT ♦ VLAMH# RHOBE » AKBE >SIGMBE»EBE * VLMDIA 
32 FORMAT ( * ANALYSIS DIDNT CONVERGE*) 

39 FORMAT ( * ANALYSIS r NO ROOT WITHIN BOUNDS' ) 

ITER=1 

RER=.001 

AERs.OOOl 

EPS=.0Q01 

IF ( NCHECK. NE. 2 ) GO TO 99 
A=10. 

B-500000 » 

GO TO 100 

99 A=.01 
B=100. 

100 IF(NCHECK.NE.l)GO TO 104 
FA=RSAF2(A> 

FB=RSAF2<B> 

GO TO 109 

104 IF ( NCHECK, NE • 0 ) GO TO 105 
FA=RSAF1(A) 

FB=RSAF1(B) 

GO TO 109 

105 FA=RSAFJU,T) 

FB=RSAF3<B»T) 

109 IF (FB*FA.GT . 0, ) GO TO 30 
GO TO 60 

38 PRINT 39 
ITER=0 
GO TO 85 

45 IF( ABS(X-B) •LToTOL). X =B* SI6N(T0L»C-6) 

A=B 

FA=FB 

IFUX~MEAN)*<B-X).GT.O.) GO TO 50 
B=MEAN 
GO TO 55 
50 B=X 

55 IFlNCHECK.NE.DGO TO 106 
FB=RSAF2(B> 

GO TO 110 

106 IF ( NCHECK. NE.O) GO TO 107 
FB-RSAF1 (B) 

GO TO 110 

107 FB=RSAF3(B#T) 

110 1TER=ITER+1 
IF(FC*FB.LT.O.) GO TO 65 

60 C=A 

FC=FA 

65 IF(ABS(FB).LT.ABS<FC) ) GO TO 70 
A=B 
FA=FB 
. B=C 
FB=FC 
C=A 
FC=FA 

70 IF(ABS(FB) .LT.EPS) 60 TO 80 
MEAN =<B*C)/2. 

IF(FB.NE.FA) X = <A*FB-B*FA)/(FB-FA) 

IFtFB.EO.FA) X = MEAN 
TOL = ABS(B*RER)4AER 
IF(ABS(MEAN-B) .LT.TOL) GO TO 80 
IF(ITER.LT.MAXITR) 60 TO 45 
PRINT 32 
ITER=-MAXITR 
80 X=B 
F=FB 

B5 RETURN 
END 
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FUNCTION RSAF3(EI*T> 

REAL L 

COMMON/ JSPl/EIEFF»FREQ»DlABOM»THlBOM»ROBOOM»WIOCE»DEGOP»WPANELrMAX 
1 1 TR » N0AT2 » ROAR AY * NW . NL r L » W * FREQR ♦ RBOOM » EBOOM » RPT » POWER ' AREA * TLEAOR 
2»BKTW»BKTL»NTYPE 

C OMMON/ JSP2/GED I AB » 6EW I DC » GEW » GET , SEE I BE * GEROBE 
COMMON/ JSP7/PCTBKL 

RL1=SQRT(T .1415S*, 2*EI*PCTBKL/(200.*T)) 

RL2=AREA/(W-( (2 . * (EI*RPT/ ( EB00M*EIEFF*4. *3. 14159) ) ** • 25 ) *12 . ) *GEWI 
1DC/6EDIAB)+TLEADR 
RSAF3=RL1-RL2 
RETURN 
END 


FUNCTION RSAFl(U) 
CCMH0N/JSP5/C0NS7 #C OEF F . EXP 
IF(W.LT .0.0) GO TO 1 
RSAF1=W-C0£FF*U** EXP-CONST 
GO TO 2 

1 R SA FI- 50 . * U 

2 RETURN 
END 


FUNCTION RSAF2IW) 

COMMON/ JSP2/GEDIAB tGEUlDC* GEW (GET .6EEIBE. GEROBE 
C OH MON /JSP6/VL AHA C> VLA MT » VLA KH tRHOBE* A KBE *SIGMBE .EBE.yLMDI A 
IF( W.LT .0.0) GO TO 1 

RSAF2-W— CW *VLAMAC *. 5-*V LA MT • .5) *• (2./3.)«( IV LAHH • ( GEW -GE ylDt )*VLMDI 
1 A* .5) / GEW ) *» (5./3» )*RHOBE/ (SORT ( AKBE *SIGMBE ) »EBE »*( 1 . /6. ) I 
GO TO 2 

1 RSAF2=50-*W 

2 RETURN 
END 


SUBROUTINE GENIGN < A.B. IGNV AL. IGNVECtN, Nu MVT R. OET B ) 
INTEGER N* NUMVTR iN 

REAL A(l), 8(1 J , IGNV AL (11 t IGNV EC ( 1 ) , DET B 
DOUBLE PRECISION DE T 
DET 8=0 . 

M=IA6S( N) 

CALL SYMDET<DET.6.8.M 
IF ( DE T . EG • 0. DC ) GOTO 1 

I F ( 0 A 8S ( DET) . GT .1 .038) D ET B=DS IGN (1 .D3 8. DET ) 

IF (DABS(DFT) .LT.l.D-38) DETB=DSIGNI 1 .D-3 8 .DE T ) 

IF( DETB.EQ.O. ) DETB=DET 
GOTO 5 

1 PRINT 3 

3 FORMA T( 25H THUS GENIGN IS UNDEFINED) 

RETURN 

5 CALL UINVRSJB.B »H ) 

CALL UTRAU(B. A, A,M»IGNVAL) 

CALL SYMIGN( A » IGNVAL . IG NVEC . N. NUMVTR ) 

IF( NUMVTR . £Q. 3 > RETURN 

CALL UXM(E,IGNVEC.IGNVEC.M»IABS( NUMVTR) 1 

RETURN 

END 


SUBROUTINE STMIGN ( At I6NV AL . IGNVEC. N. NUMVTR ) 

PARAMETER MXN=100 

INTEGER N*NUHVTR.NPi» I#.» . H » NH M, IM1 , MM1 . K, KP1 , L AST . L . F. J 
REAL A( 1) tIGNVAL(D»IGNVEC<l)»NORM*OVRFLO*EPG»LAMOA»DIAtSCALEf 
1 C. D. DO. FF.EPS.PRECS»S*T*XX»DI(HXN). C0( MX N> . CD Q( MX N) 

LOGICAL FULL 

C0HH0N/SR1/ 0I/SR2/ CD/SR3/C00/ NR M/ NORM 
CALL TRIDIA( A .01 .CD » COG «N ORM * I AB S ( N ) ) 

FULL=.TRUE. 

GOTO 10 

ENTRY TRIGN ( IGNV AL • IGNV EC » N, NU MVT R ) 

FULL=. FALSE. 

10 PRECS =1. E- 8 
OV RFLO- 1. E38 
L=IABS(NI 

SCALE=SGRT(OVRFLO*PRECS )/ AMA Xl(l. t A MINI ( NORM . OV RFLO* PRECS J) 

S-S CALE* *2 

EPS= PRECS* NORM* SC ALE /FL CA T( L ) * * 2 
EPQ = AMAX1( EPS**2. 1000. /OVR FLO) 

NP1=L*1 
00 IS 1=1. L 


15 

P=NP1-I 

IGNVALIP >=DH I) 
DI(I) =01(1)* SC ALE 
COQ(I) =CDQ(I)*S 


20 

M=L 

GOTO 2* 

LAHDA=DI( M)/SCALE 


21 

COQ(M) =IGNV AL (NNM ) 
I=NMM 

IF ( I.EQ.l ) GOTO 22 



I M1=I -1 

IF (LAMDA.GE.IGNVAL( IM1II 

GOTO 22 


IGNVALd >=IGNVAL (IM1) 

I=IM1 

60T0 21 
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22 IGNVAM U -L AHDA 

IF ( M- EQ • 1 > GOTO 99 
M -HMl 
CDQ ( M )-0. 

29 MM1=M-1 

IFIH.E0.1) GOTO 20 
25 K=H 

30 K-K- 1 

IF ( K.E Q. C ) GOTO 35 
IFtCDQIK) .GT.EPQi GOTO 30 
IF(K.EQ.HMl) GOTO 20 
CDGI K) -0. 

35 K-K^l 

T=OItH>-DIINMl > 

XX=CO Q I H Ml ) 

DD=T/2* 

IFIDO'DD.LE.ABSI PRECS*XX H GOTO 37 
LAMDA=IXX/DD)/tl.*S3RT ll-*XX/0D**2l I 
GOTO 90 

37 LAMDA=SQRTl XX) 

90 D=DII M)*LAMDA 

DD-DI ( MH1 )-LAM0A 
IFIK.NE. MMl) GOTO 92 
OIIN)=D 
OltMHl »=00 
GOTO 20 

92 LAMDA=D 

IF t A8S (T) .LT.EPS.AND .ABS(DD) . LT . A BS ( 0) ) L AMD A=0 0 
DIA=OI(K >-LAHDA 

IF ( ABS ( 01 A) -LT.EPS1 DI A-D I A^SIG N( E PS* D I A ) 

FF=DIA**2 

XX=COQ IK) 

T=XX*FF 
GOTO 99 

93 FF=DIA/C*DIA 
XX=CDQt KP1I 
T=XX*FF 
CDQl K>=S*T 
K=KP1 

99 KP1=K»1 

C-F F/T 
S=XX/T 
XXrOIA 

01 A=C*I 0KKPD-LA1 DA)-S*XX 
Dll K )= < XX-0 IA I ♦DI IK PI I 

IF| ABSI DIA) «LT «EPS ) DIA=DIA»SIGNI C*EPSt OIAI 
IF(K.LT.HN1» 6 0 TO 93 
COOI MN1 I- 01 A/ C • DI A *S 
DIIM>=OIA*LAMDA 
GOTO 25 

99 LAST=NUHVTR 

IFIN.GT.O) GOTO 55 

IfilGHVALtl) .LT.O..AND.IGNVALIL».GT.O.> 60T0 50 
IF|ISIGN|1*LAST).GT.0J GOTO 56 
LA ST= -LAST 
GOTO 57 

50 F = 1 
K-2 

IFIIGNVALI2J .GT.O.I GOTO 59 
K=3 

51 IFIIGNVALIKI .GT.O.) GOTO 52 
K=K*1 

GOTO 51 

52 J=K-1 

P" I K~F ) /2*F-1 
DO S3 I = F*P 
DD=I6MVAUI) 

IGNVALlI)=IGNVAL I J> 

IGNVALI J >=00 

53 J=J-1 

IF(K.GE.L) GOTO 55 

59 F -K 
K- NP1 
GOTO 52 

55 IF IISI6NI 1 .LAST I .GT • 0 J GOTO 57 
L AST --LAST 

56 F- 1 
K=NP1 
GOTO 52 

57 IFINUMVTR.EG-OI GOTO 70 

P=0 

00 65 1=1* L AST 

CALL TRIVECI CDS *CD *01 *IGNVAL»I»NORM»L) 

IFIFULLI CALL TR IBAK I A . C 0 » D It L» 

H=P 

DO 60 K=1*L 
H-M+l 

60 I6NVCCIH) =DIIK I 

65 P-P* L 

70 IFIN.GT.OI RETURN 

DO 75 1=1 tL 

75 IGNVAL(II=l ./IGNVAL«I» 

RETURN 

END 
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SUBROUTINE MBIM.W.L . M. N2.N3.N4.N) 

INTEGER NX » N2 . N 3 » N4 »N 

REAL KIN.ll.W.L 

A=M/32. 2/420- 

B- A*L 

C=B*L 

Ml N1.M1 )=HIN1.N1 )*156.*A 
MIN1.N2I =M<Ni.M2»*22.*8 
K(N1*N31=K(N1.N31*S4-*A 
MINI* N 4) =H I N1 *N4I-13.*B 
Ml N2 »N1 I=KIN2*N1 )+22. *B 
MCN2»N2I "M I N2 *N2)+4.*C 
Ml N2.N3 ):HCN2.N3 )*13.*B 
M(N2.N4)=HlN2»N4)-3.*C 
Ml N3.N1 )=MI N3» N1 )+54 • *A 
MIN3»N2)-M(N3»N2) + 13-*B 
Ml N3.N3 >rKCN3»N3 1*156 - *A 
MIN3.N4I =HIN3*N4|-22.*B 
Kl N4 iNl )=>IN4»Nl )-13. *B 
MIN4.N2) =N<N4 iN2)-3.*C 
HI N4 .N3 )=KCN4.M3 >-22.*B 
M(N4.N4I=M(N4»N4)*4.*C 
RETURN 
END 


SUBROUTINE SYHDET lOET.A.UtNI 

INTEGER N> CtCCtLCCt CPI $ COL tRCtR COL* DIAGt CCPlt 0IAGM1 
REAL A I 1 ) »Ut 1) .SUM 
DOUBLE PRECISION OET 
SUM=AI1 ) 

IF| SUN-LE- 0- I GOTO 10 
U« 1I=SQRTISUM) 

DET=UI1I 
GOTO 12 
10 PRINT 1 

1 FORHATI ////4 7H SYMOET UNDEFINEO* MATRIX NOT POSITIVE DEFINITE I 

IF< SUM.EQ.O.) PRINT 2 

2 F0RHATI1H**46X*17H. DE TER MIN ANT =0-0) 

OET =0.00 

RETURN 
12 RC=2 

DO IS C=2*N 

UIRC) = AI RO/Ulll 
15 RC=RC*C 

CC=1 

DO 35 C=2.N 
QI AG=CC* C 
CCPl=CC*l 
0 1 AGM1 =0 1 AG- 1 
SUM=0. 

DO 20 RC=CCP1 * DI AGH1 
20 SUM=SUH* U( RC ) **2 

SUM=AI DIAGI-SUM 
IFISUK.LE.O. I GOTO 10 
Ul DIAG|=SORTI SUM ) 

DET=DBLEIUCDIAG)I*DET 
IFIC.EQ.N) GOTO 35 
CP1=C*1 
LCC=0IAG»1 
DO 30 C0L=CP1*N 
SUM=0. 

RCOL= LCC 

00 25 RC=CCP1. DIA6M1 
SUM=SUH*U< RCOL > *U IRC > 

25 RC0L=RC0L*1 

Ul RCOL )= (A I RC OL ) -SUM) /U I DIAG I 
30 LCC=LCC+ COL 

35 CC=DIAG 

IFIOABSC DETI.LT.1.D-153J GOTO 40 
IFIDABSIDET) .GT. 1.0*153 > GOTO 45 
DET=DET**2 
RETURN 

40 WRITE 1 6 *41 1 

41 FORHATI//* SYMOET UNDERFLOW* DETERMINATE SET TO l.D-306 *) 
DtT=l.D-306 

RETURN 

45' WRITE 1 6 *46 ) 

46 FoRMATI//' SYMOET OVERFLOW* DETERMINATE SET TO 1.0 *306M 

0ET=1.0306 
RETURN 
ENO 


SUBROUTINE SYMINV IA.AI.N.DET) 

INTEGER NfJ*N* NM1.R .RP1.R I. I* IC. C. CM1. DIAG. ROW . ROW I* ROW R.FLAG 
REAL Af NI.AKNI.SUM.OET 
OOUBLC PRECISION DTR 
CALL SYMDETI DTR.A.AI.N) 

FLAG=1 
or t=o • 

IF(DTR-EQ.O.OO) GOTO 1 

IF (DABSCOTRl .GT.1.036 I OET =DSIGN( 1.D38 . 0 TR 1 
IFIOABSI DTRI.LT. 1.0-38) DET = OS IGN 1 1. D- 38. DTR > 

IF IDE T.EG-O. ) OET =OTR 
GOTO 5 

1 PRINT 3 

3 FORM AT I 2 5H THUS SYMINV IS UNDEFINED) 
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RF TURK 

ENTRY SIMVRS (AI.NJ 

FLAG=1 

GOTO 5 

ENTRY UINVRS IA.AI.N) 
feag=o 

S KN=IN**2*N)/2 
DIAG=0 
KM1=N-1 
DC 3D R=l. NHL 
DTAG=DIAG*R 
AI ( DI A G I -1. / All DIAGJ 
PP1 = R-»1 
IC=DIAC*R 
DO 20 C=RPliN 
RI=DIAG 
SUH-Q . 

CM1=C-1 
CO 1C I-R.CK1 

5UM=SUK *AI(Rl)*AII IC ) 
RI-RI* I 
10 IC=IC*1 

All RI >= -SUM/A II I C ) 

IC=RI*C 

20 CONTINUE 

30 CONTINUE 

AI CNN|=1./AI (MN) 
lF|fLAG.E3*0> RETURN 
DI AG- G 
DO GO R=i*N 
ROUR=DIAG-H 
D I AG -DI AG.R 
DC 50 R0W=1.R 
SUM=0. 

RI-DIAG 
ROUI=ROWR 
CO 40 I-R.S 

SU K-S UM ♦ All ROW I > .A IlRl> 
ROUI- RCUI+ I 


43 

RI=RX»I 


All ROWR > = : 

50 

R0UR=R0UR»1 

GO 

CONTINUE 


RETURN 


CMC 


SUBROUTINE K RA VI K »N »E IB EA M t E IB OOH » TOK * TI K *T OB »T IB . NW • NL » L »W * SYM > 
INTEGER N» NL»NN.NW .NZ»NW C# NBN, I. JtNl. SYM 

REAL KC II »EI »E IBOOK t EIB EA M t T OM * TIN »T OB t T IB iOTB *0TM tL * LL • LU »U • T 
EI=EI800M/2. 

NUC- N W+ 1 

NZ =NW C*NL* I 2* NW C) 

NB N- NWC • I NL«1) 

LW = U/FU0AT (NW ) 

LL=L/FL0ATINL1 

otm - i tom-timi /Float inl w> 

DTB- 1 TI6-TCB >/FL OAT I NL ) /2 . 

T=TIM/FL0ATINW) ♦DTM/2. 

NN=NBN* 1 
N =N Z • NZ 
DC 5 1=1. K 
5 KIU-G. 

DC 10 1 = 1 .NW 

10 CALL KBIK.EIBEA-.LW.NM-I.I. NBN- I. I+l.NZJ 

DO 15 1 = 2. NWC 

15 CALL KSKlK.LL .T.MZ.NU* I. NUC^I.NZ. NZ) 

K1=NUC+ 1 
NN=NBN-NWC 
DC 20 J= N1 » N N » NWC 
. T =T + DT M 
DO 20 1=1. NW 

20 CALL KSMIK.LL.T* J*I- 1. J«-J +NW . J* I+NWC. J>I. NZ > 

IF I SYM. EG* 1 ) GOTO 30 

Nl=NBN+3 

T- I TOB-DTB 1/2. 

DC 25 I=N1.NZ .2 
T=T*DTE 

CALL KSBIk. T.LL. I- 1 . I. I-3.I-2.NZ> 

25 CALL KBIK.EI .UL.I-l.I .1-3. I-2.NZ1 

CALL KSBIK . T* DT B.LL » N Z » N Z. NZ- 2 . NZ- 1 . NZ ) 

CA LL KBIK.EI.LL.NZ.NZ.N Z- 2 . N Z-l » NZ ) 

CALL PAKSJeiK .K .2.NZ-1 .MZ> 

N= NZ-2 
RETURN 

30 CALL PAKSUBlK.K.l.NBN-1 »NZ> 

N=NBN-1 

RETURN 

END 
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SUBROUTINE TR IDIA ( At D IA , CO DI » COSO* NORM* N) 

INTEGER NN«N «NH1 *ItIHl*COLtC tC P 1 »CCP 1. LCC • ROW * R* 

1 RP1» RP2 tRC.OIAG. DI AGPR 

REAL At N) .CODKN ) .CDSG(N) »DIA t N ) .NORM»SUM.W t X 
NORM -0- 
NN=C N**2+N)/2 
NM1 =N-1 
CCP1-1 

DO 20 C-l » NM1 
CF1=C*1 
DIAG=CCP1*C-1 
SUM- 0 . 

DO 10 RCrCCPltOIAG 
10 SUN=SUH*AES< A < RC > > 

RC=DIAG*C 
DO 15 C0L=CP1.N 

SUH=SUM+ABSI A < RC ) ) 

15 RC=RC*COL 

IF(SUN.GT.NORM) N0RM=5UH 
CCP1=CCP1*C 
20 CONTINUE 

SUH=0. 

DC 25 RC=CCPl»N'N 
25 SUM=SUH*AE St A t RC ) J 

IFlSUM.GT.NORM) NORM=SUl 
DIAGEO. 

3D 00 95 R-l » NM1 

RP 1 - R*1 
DIAG=DIAG*R 
DIA ( R A (D IAG ) 

DIAGPR=DIAG*R 

SUK=0. 

RC=DIAGPR 
DO 35 C=RP1.N 

SUH-SUM*AtRC)**2 
35 RC- RC*C 

CDSQ<R»=SUM 
CODI( RI^SORT (SUM) 

IF( NORM- CODICRJ.NE. NOR M) GOTO 40 
COOI < R) = G. 

CDSQCR)=0. 

GOTO 95 

40 I F ( A( 0 1 AGPR ) . GT . 0 . ) Cl Cl <R ) =-CODI tR) 

A(DTAGPR) = A(DIAGPR)-CCDIIR ) 

W - At DI AGPR) *CODItR) 

RC = DIAGPR 
DO 50 C-RP1 »N 
DIA 1C )rA (RC ) 

SO RC=RC*C 

LCC-DIA GPR *1 

RC=LCC 

SUM=0. 

00 55 C-RP1 >N 

SUM-SUM-* A t RC ) *D I A <C ) 

55 RC=RC*C 

COOK RP1) = SUM/W 
RP2=R*2 
LCC- LCC * RP 1 
IFtR.EQ.NMl) GOTO 71 
DO 70 I=RP2tN 
IM1-I- 1 
RC- LC C 
SUM-0. 

CO 60 ROUrRPltlMl 

SUM-SUM*AtRC) * D I A( RO W ) 

60 RC=RC*1 

DO 65 C=I»N 

SUK-SUK*A( RC ) *DIA CC ) 

65 RC=RC*C 

CCDI(I)=SUM/W 

LCC=LCC*I 

70 CONTINUE 

71 SUM-0. 

DO 75 I=RPltN 

75 SUM = SUM*CIA t I) * CODIt I) 

X-SUM/W/2. 

DO 80 I-RP1 »N 

80 CODHI ) = CIA( I )*X*CODI (I ) 

LCC=OI AGPR* 1 
DO 9C C — RP1 1 N 
RC=LCC 

DO 85 ROU= RP1»C 

A( RC) =DIA(ROW)* COO It C ) * CODI < R OU ) * DI A t C) ♦ At RC) 
85 RC- RC ♦ 1 

LCC=LCC*C 
90 CONTINUE 

95 CONTINUE 

COCK N )"C . 

COS Q ( N ) - 0 . 

D I A ( N ) = A ( NN ) 

RETURN 

END 
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SUBROUTINE TRIVECIDI.COOI.V.IGNVAL. NuM.NORh.N) 

INTEGER N»NUM,ItIPl*J»NMl »ITR 

REAL OZ(N> .COOICN) »V(N)t IGNVALIN) . X .X 1 . X2 • Y »Z * SU Hr C, E. S# 
1 PI IOC) #OI1CO).WI1CO )» NORM t EPS. PRECS 

LOGICAL INTNHI » InT 1100 ) 

C0KMGN/SR4/P/SR5/G/SR6/INT/SR7/U 

PRECS =1 • E- 8 

Xl-D. 

X2=S3RTI 2. ) 

NK1=N-1 

E=IGNVAHNUH> 

IF t NUN. EG. 1 ) GOTO 10 

IF( ABS| E-IGNV Al (K'JM-1) ). LT -E.pRECS* 100. > GOTO 15 
10 DO 12 1 = 1. N 

12 V(I)=1. 

GOTO 25 

15 DO 20 1=1. N 

X= A MOD I XI* X2 .2. ) 

X 1 =X2 
X2= X 

20 V( II =X-1 . 

25 E PS= PRECS *NC RM/FLOA T I N ) •• 2 

Y=l. 

X= IDI 1 1 )-E l/NORH 
DO 35 1=1. NH1 
IP1=I*1 

S=CODI(T)/NORM 
C=IDI IIP1 ) -E ) /NORM 
Y =S • Y 

IFIABS(X).GE.ABSIS) ) CO TO 30 
P 1 1 ) =S 
Q I I )=C 

INTI I) r.TRU E. 

Z=-X/S 
X=Y*Z* C 
Y= Z 

GOTO 35 

30 IFIABSI XI.LT.EPS) X=E FS 

P I I ) =X 
GII )= Y 

INTIIJ =• FALSE* 

Z=-S/X 
X =C*Z* Y 
Y=l. 

35 WII)=Z 

INTNK1=INTINK1 ) 

IFIABSIXJ.LT.EPS) X=EPS 
I T R= 0 

40 ITR=IT R. 1 

VIN)=VIN)/X 
INTINM1) =. FALSE. 

SUK= V I N ) * • 2 
DO 45 J=1.NM1 
1 = N-J 

Y=VI I ) - G C I)*VI 1*1) 

IF I INT 1 1 I) Y = Y-C ODI 11*1) • V 1 1*2 > / NCRM 
VI I) =Y/P( I) 

45 SUK=SUM-» VI I) **2 

S =S 3RT (SUM) 

DO 50 1 = 1. N 
50 V(I)=VII)/S 

IFIITR.GE.2) RETURN 
INTI NH1) =INTNK1 
DO 60 1=1. NK1 
IP1=I*1 
Z=VII ) 

IF ( INT ( I) ) GOTO 55 
VI IP1 ) = V IIP1 )*Z* VII ) 

GOTO 60 

55 Y=VII) 

VI I) =V I IP1) 

VIIP1 )= Y *Z • V 1 1 ) 

60 CONTINUE 

GOTO 40 
END 


SUBROUTINE CP YSJBlH.SU 8. Ri »RZ# Cl »C2.N) 

INTEGER Rl .R2.C1 .C2.R1C 1.R1C2.R1C .R2C »RC .IJ.N.NN 
REAL Nil I.SU3I1) 

R1C1 = N* 1C 1-1 )+ Rl 
R1C2 =N • C2 
R2C=R1C1+R2-R1 
NN = N 
GOTO 10 

ENTRY PAKSU8IH.SUB.R1.R2.N » 

R1C 1 = N* I P.1—1 )*RI 
R1C2=N*R2 
R?C= R1C 1 
NN=N* 1 
10 IJ=0 

DO 20 R1C=R1C1 .R1C2.N 
DO 15 RC=R 1C . R2C 
IJ=IJ+1 

15 SUB 1 1 J )= HI RC ) 

20 R 2 C=R2 C* NN 

RETURN 
ENO 
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SUBROUTINE KRA Y I H »N » ROB EA M»R 00 OGM » ROA RA Y . NW . NL . L • U .S Y K ) 
INTEGER N. NL. NN.NU .NZ.NW C, NBN » I* J.N1 , SYM 
REAL Ml 1) .ROEEAK.ROBCOM »R OAR A Y t WO . WE . UA • L tLU . LL . U 
MWC-NW + 1 

NZ = NUC-»NL*I2+NWC ) 

N BN =NWC* I NL*1 ) 

LW-W/FLCA Tl NW ) 

LL = L/FLO AT I NL ) 

WO=ROSOOM*LL/2. 

WF=ROBEAM*LU 
UA = RCARAY*Ll *LW 
NN=N8N+1 
N=NZ*NZ 
00 5 I— I * N 

5 Mil ) - 0. 

DC 1C 1=1. NW 

10 CALL MBIM.WE .LW.NN-I. I. NBN -I • I + l »NZ » 

00 15 1=2. NWC 

15 CALL MKEHI H.WA .NZ.NW+I. NWC -*1 , NZ . KZ J 

N1 =NWC+ 1 
NN = NbN-NWC 
00 20 J=N1 » NN .NWC 
DO 20 1=1 . NU 

20 CALL MM EM (K.WA.J^I-l » J* I* NW . J+I +NU C. J+ I . NZ ) 

iriSYM .EQ. 1) GOTO 3C 
Nl=NBN«-3 
DO 25 I=N1.NZ*2 

25 CALL M8l M.WO.LL » I“l» I. I-3.I-2.NZ) 

CALL MB ( H .WO »LL . NZ.N Z .N Z- 2 . N Z-l . NZ > 

CALL PAjKSUBlM .Mf2.NZ-l.NZ) 

N= NZ-2 
RETURN 

30 CALL PAKSU8 I M.H.l.NBN-1 »NZ> 

N=N BN- 1 
RE TURN 
END 


SUBROUTINE UTRAU IU. A . J T AU » N . SR AT CH) 

INTEGER N.CtCC tCCPR.LCC »C OL . R . RC » RCPI . RI » I • I PI . DI AG 
REAL AINJfUiN I .UTAU(N) »S RAT CHIN) .SUM 
CC = IN**2-»N)/2 


10 CC=CC-C 
DI AG=0 
DO 40 1 = 1. C 
LCC= DI AG 
DIAG=DIAG-*I 
RI=L CC 
RC = CC 
SUM = C. 

20 RZ=RI + 1 

SUM = SUK*A( R I ) * U I R C ) 
IFIRI.LT.OIAG) GOTO 20 
IFII.EG.C) GOTO 35 

IP1 = I*1 
RI = RI*I 

00 3C COL =IP1 » C 
RC = RC ♦ 1 

S'JK=SUM + A(RI) * U ( RC ) 

30 RI=RI*COL 

35 SRATCHII) =SUM 

40 CONTINUE 

RC=CC*C 
R=C 

45 RC =R C- R 

SUM=0. 

DO 5G 1=1. R 
RCPT=RC+I 

50 SUM=S UM *U IRCPI) • SR AT CHI I) 

CCPR=CC ♦ R 
UTAUICCPR )=SUM 
R= R -1 

IFIR.CT.O) GOTO 45 
C = C-1 

IFIC.GT.l) GOTO 10 
UTAUI 1 )=AI1) *U 1 1 ) • *2 
RETURN 
END 


SUBROUTINE U XM I U . M » UM . N »P ) 

INTEGER N.P.NP.NMl.NC.OC.CC.R CC.R I. IC 
REAL Ul N) .Ml N) .UMIN) .SUM 
NH1 =N-1 
NP = N*P 

DO 30 0C=1.NP.N 
NC = CC ♦ MMl 
R C C=1 

DO 20 C C = CC • NC 
RCC=RCC+CC-OC 
RI= RCC 
SUM=0. 

DO 10 IC =C C » NC 
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RI=RI+IC-OC 

ID SUK-SU K*UIRI)*K(IC) 

UM|CC» -SUM 

2D continue: 

30 CONTINUE 
RE TURN 
END 


SUBROUTINE TRIBAK IA.CODI.VtN> 
INTEGER Nf NH1 t C* I.R.RP1. RC.DIAGPR 
REAL A ( N ) .COC-I I N ) .V I N ) » SUM.F *X 
NM1 =N-1 

GIAGPR=IN**Z*NJ/2-l 

R=N 

DO 23 1=1 *t«Hl 
R=N-I 
RC=DIAGFR 
RP1=R+1 

IFICODI IRJ.EQ.C.) GOT C 2D 
X=AIRC> 

SUM=0. 

DC 10 C=RP1 » N 

SUH-SUK* A < RC I ♦ V I C 1 
10 RC=R C* C 

F=SUH/< X*C CD I ( R I ) 

RC-OI AGPR 
DO 15 C = RP a » N 

V<C>=AIRC>*F+VIC) 

15 RC=RC*C 

20 DIAGpR = DIAGPR-RP1 

RETURN 
END 


SUBROUTINE MM EH (K * W.Nl »V 2# N3. N4. N) 

INTEGER N1.N2.N3.N4 .N 

REAL MIN.ll.W 

A= U/3 6 . /? 2. 2 

8=2. ♦ A 

C- 2. *B 

MINl.Nl) =MlNl *N1) +C 
Ml N1 »N2 )= Ml Nit N2 MB 
MINI »N3) =M(N1 »N3> ♦ A 
Hf Nit >=KI Nit *ll» )*8 
K(N2»N1)=M<N2.N1>*B 
Ml N2tN2 1= Ml N2 » N2 >*C 
MIN2.N3) =H I N2 . N 3 ) ♦ 8 
M l N2 »N4 )=NIN2»N4 HA 
N I N3* N1 ) =H I N3 » N1 ) ♦ A 
Ml N3 t N2 >=M(N3»N2 )*B 
MIN3»N3) =M IN3.N3J+C 
Ml N3 t N4 )= Ml N3. N4 )*B 
MIN4»N1)-MIN4»N1) + B 
Ml N4.N2 )= M I N4 t N2 MA 
M (N4 »N3) =M(N4 .N3)*B 
Ml N4tN4 )=KIN4.N4 MC 
RETURN 
END 


SUBROUTINE KEIK.EI.L.N1 .NZ.N3.N4.N) 

INTEGER N1 tN2 «N3tN4tN 

REAL KINtDtEItL 

A=2. * EI/L 

B=3.*A/L 

C=2. * B/L 

Kl N1.N1 1=K(N1»N1 )*C 
K IN 1 » M2 ) =K ( N1 t N2 ) ♦ B 
Kl NltN3 )=KINltN3 )-C 
K INI »N4) =K I NL »N4)*B 
Kl N2 rNl )=KI N2.Nl ) + B 
K I N 2 • N 2 ) =K I N2 tN2)+2.*A 
Kl M2.M3 J=KIN2»N3 1*6 
KIN2.N4 J =K (N2 »N4) ♦ A 
K I N3 t N1 1= K ( N3 t N1 I -C 
K IN3.N2) =K f N3 tN2 ) - B 
K I N3 • K3 J- Kl N3» N3 l+C 
K IN3tN4) =K I N3 » N4) - B 
KIN4.Nl J=K( N4.N1 )♦& 

KIN4.N2) =K IN4 »N2 > ♦ A 

Kl N4 »N3 )=K I N4 »N3 ) -B 

K I N 4 . N 4 ) =K (N4 » N 4 ) *2.* A 

RETURN 

END 
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SUBROUTINE KSHCK»L »T»N1»N2»N3»N4»N) 
INTEGER N1»N2»N3*N4»N 
REAL T»L*K(Ntl) 

A - T/6./L 
B -2 • * A 

Kt NI » NI >*e 

K ( N 1 • N 2 ) -K ( N1 » N 2 ) - B 
K t NI t N3 )- K ( N1 » N3 ) -A 
K(Nl»N<l)=KtNl»N4)«A 
K I N2 » N1 ) — K < N2t NI ) -B 
K(N2*N2)=K(N2»N2)+B 
KtN2»N3 »=K(N2*N3 t*A 
KIN2»N4 ) =K (N2 *N4) - A 
M N3 » Ml ):K( N3 v N1 ) -A 
K t N 3 » N2 t -K t N3 t N2 ) *A 
Kl N3t N3 >=K<N3.N3 1 *B 
KtN3».N4)=K (N3»N4)-B 
Kt N4 » NI JrKfN^tNX >*A - 

KIN < 4tN2> =K<N4 » N 2 ) - A 
Kt N4 ,N3 )=KI N4tN3 ) -6 
K ( N 4 » N 4 ) =K (N4 »N4J +B 
RETURN 
END 


SUBROUTINE KS B(K#T »L»N l.N2f N3.N4. N) 

INTEGER NI * N2» N3 »N4 »N 

REAL K < N • 1 ) *T *L 

A3 6* • T/5. /L 

B3T/1D. 

C= T*L/3D. 

K(N1»N1)3K(NI»N1)+A 
Kt NI *N2 )3 K t N 1 1 N2 )+6 
K(NltN3l=K(Nl*N3)-A 
Kt N1.N4 I3K1N1.N4 ) *B 
K(N2iNll 3K(N2.Nlt ♦B 
Kt N2*N2 I3KtN2»N2 )*4. *C 
K ( N2 i N3 ) 3K1N2 tN3>-B 
Kl N2.N4 >= K ( N2 t N4 ) -C 
KtN3»Nl)=KfN3»Nl)-A 
Kt N3.N2 )=KtN3.N2 ) -3 
KtN3»N3> 3KtN3*N3) +A 
Kt N3 » N4 )=Kt N3.N4 ) -B 
K(N4.Nl)3KtN4,Nl)*B 
KIN4.M2 >3 K t N4 *N2 >-C 
KtN4.N3)=KtN4 .N3J-E 
Kt N4 # N'4 J3K(N4*N4 >*4. *C 
RETURN 
END 
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